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ABSTRACT 


Steam  condensation  heat  transfer  measurements  were  made 
in  a  5-tube  test  condenser  having  an  additional  perforated 
tube  to  simulate  up  to  30  active  tubes.  Results  were 
obtained  for  smooth  tubes  and  roped  tubes  wrapped  with  wire 
A  Sieder-Tate  equation  was  used  to  correlate  the  inside 
heat-transfer  coefficient.  For  smooth  tubes,  a  leading 
coefficient  of  0.029  was  found,  while  it  was  0.061  for  the 
roped  tubes.  The  average  condensing  coefficient  measured 
for  30  smooth  tubes  was  0.59  times  the  Nusselt  coefficient 
calculated  for  the  first  tube.  When  the  smooth  tubes  were 
wrapped  with  wire,  this  ratio  increased  up  to  0.86. 

Further,  roped  tubes  without  wire  experienced  a  ratio  of 
0.63,  while  roped  tubes  wrapped  with  wire  resulted  in  a 
ratio  of  0.86.  These  preliminary  data  show  that  wire- 
wrapped  tubes  may  lead  to  a  significant  reduction  in 
condenser  surface  area. 
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NOMENCLATURE 


2 

Aq  Outside,  heat-transf er  area  of  one  tube  (m  ) 

2 

A^  Inside,  heat-transfer  area  of  one  tube  (m  ) 
Sieder-Tate  coefficient 

CpW  Specific  heat  of  water  evaluated  at  (KJ/kg*K) 

0  14 

Cf  Correction  factor  ( n/|aw ) 

Inner  diameter  of  the  tube  (m) 

D  Outer  diameter  of  the  tube  (m) 

o 

2 

g  Acceleration  of  gravity  (9.81  m/s  ) 

h.  Experimentally-determined  value  for  the  inside,  heat- 
transfer  coefficient  (W/m2*K) 

hf  Latent  heat  of  vaporization  (KJ/kg) 

h  Local,  outside,  heat-transfer  coefficient  for  the  Nth 
tube  (W/m2K) 

h  Heat-transfer  coefficient  calculated  from  the  Nusselt 

u  equation  (W/m2*K) 

h,  Outside  heat-transfer  coefficient  for  the  first  tube 
( W/m2 • K ) 

k^  Thermal  conductivity  of  the  condensate  film  (W/m*K) 

k  Thermal  condictivity  of  the  cooling  water  evaluated 

at  Tfa  ( W/m* K ) 

km  Thermal  conductivity  of  titanium  (W/m.R) 

L  Condensing  length  (m) 

LMTD  Logarithmic  Mean  Temperature  Difference  (JC) 
m  Slope  of  the  least-squares-f it ,  straight  line 

m  Mass  flow  rate  of  cooling  water  (kg/min) 
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The  number  of  tubes  in  a  column  or  the  tube  number  of 
a  given  tube 

Water-side  Nusselt  number 
Prandtl  number  evaluated  at 
Heat  transfer  rate  (W) 

2 

Heat  flux  based  on  outside  area  (W/m  ) 

Water-side  Reynolds  number 

2 

Fouling  thermal  resistance  (m  K/W) 

Wall  thermal  resistance  based  on  the  outside  area 
(m2K/W) 

Outside,  local,  heat-tranfer  coefficient  ratio  (i^/h^) 
Outside,  average,  heat-tranfer  coefficient  ratio 

<Vhi> 

Spacing-to-diameter  ratio  of  tubes 
Average  cooling  water  bulk  temperature  (°C) 

Cooling  water  inlet  temperature  (°C) 

Cooling  water  outlet  temperature  (°C) 

Wall  temperature  (°C) 

Average  condensate  film  temperature  (°C) 

Saturation  temperature  of  steam  (°C) 

Vapor  (steam)  temperature  (°C) 

2 

Overall  heat-transfer  coefficient  (m  K/W) 

.  2 
Outside  heat-transfer  coefficient  (m  K/W) 

Cooling  water  velocity  (m/s) 

Sieder-Tate  parameter  (X  =  Re+0  *  ^Pr^”^^^)  (— )  ^  * 

yw 

Temperature  difference  (T  -T,  )  (°C) 


5  Heat  capacity  parameter  (C  AT/h,  ) 

Er  J 

2 

H  Dynamic  viscosity  of  water  evaluated  at  (N*s/m  ) 

2 

Dynamic  viscosity  of  condensate  evaluated  at  (N*s/m  ) 

2 

u  Dynamic  viscosity  of  water  evaluated  at  T  (N-s/m  ) 

w  w 

C  Acceleration  parameter  (KAT/ph^g) 

p  Density  of  cooling  water  (kg/m^) 

Pv  Vapor  density  (kg/m^) 
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In  recent  years,  there  has  been  a  continued  interest  in 
the  reduction  of  the  size  and  weight  of  propulsion  systems 
aboard  both  surface  vessels  and  submarines.  Especially, 
actual  dimensions  of  naval  condensers  have  a  critical  bear¬ 
ing  on  cost  and  performance  of  the  ship.  Often,  compactness 
is  more  important  than  thermal  effectiveness  when  the  over¬ 
all  performance  of  the  ship  is  considered;  in  a  submarine, 
the  diameter  of  the  pressure  hull  can  depend  on  the  dimen¬ 
sions  of  the  condenser. 

The  importance  of  compactness  justifies  measures  to 
raise  the  overall  heat-transfer  coefficient  of  condenser 
tubes  despite  the  penalties  which  may  occur,  i.e.,  the 
increased  pumping  power  and  tube  cost. 

Naval  condenser  design  is  based  upon  the  Heat  Exchange 
Institute  (HEI)  specifications  for  steam  condensers  [Ref.  1] 
and  also  the  standards  of  the  Tubular  Exchange  Manufacturers 
Association  (TEMA)  [Ref.  2],  Search  [Ref.  3]  investigated 
the  present  condenser  design  processes,  including  the 
feasibility  of  enhanced  heat  transfer  in  naval  condensers. 

He  concluded  that  the  current  design  is  very  conservative, 
and  he  predicted  that  a  forty-percent  reduction  in  condenser 
weight  and  volume  could  be  achieved  depending  on  the  heat- 
transfer  enhancement  method  used. 


In  recent  years,  many  research  efforts  have  been 
directed  to  the  study  of  heat-transfer  enhancement  tech¬ 
niques  and  their  application  to  heat-exchanger  design.  Webb 
[Ref.  4]  has  summarized  extensive  works  of  augmentation 
techniques.  At  the  Naval  Postgraduate  School,  Beck  [Ref. 

5],  Pence  [Ref.  6],  Reilly  [Ref.  7],  Fenner  [Ref.  8]  and 
Ciftci  [Ref.  9]  conducted  experimental  research  into  various 
kinds  of  enhancement  schemes  employing  a  single-tube  test 
condenser. 

The  above-mentioned  investigations  concluded  that,  for 
the  same  diameter  tube,  the  overall  heat-transfer  coeffi¬ 
cient  of  enhanced  tubes  can  exceed  those  for  smooth  tubes  by 
almost  100  percent.  Reilly  and  Fenner  [Ref.  10]  revealed 
that  most  of  the  above  mentioned  augmentation  occurred  on 
the  cooling-water  side  due  to  a  combination  of  increased 
surface  area,  and  increased  turbulence  and  swirl  in  the 
cooling  water  flow.  Little  or  no  improvement  occurred  on 
the  steam  side.  Eissenberg  [Ref.  11]  performed  an  extensive 
study  on  condense r- tube ,  heat-transfer  coefficients  using  a 
multi-tube  bundle. 

In  order  to  investigate  the  outside  heat-transfer 
performance  of  various  enhanced  tubes  in  tube  bundles, 
research  was  conducted  at  the  Naval  Postgraduate  School. 
Noftz  [Ref.  12]  modified  a  test  apparatus  initially  designed 
by  Morrison  [Ref.  13]  to  simulate  a  tube  bundle  using  five 


active  tubes  arranged  in  a  vertical  plane.  A  perforated 
tube  was  located  at  the  top  of  the  bundle,  through  which 
water  was  flooded  to  simulate  bundles  having  up  to  30  tubes 
in  a  vertical  row.  His  investigations  determined  that  the 
heat-transfer  coefficients  for  a  given  tube  of  the  tube 
bundle  increased  as  the  mean  vapor  velocity  increased,  but 
decreased  as  the  amount  of  condensate  inundation  increased. 
The  experimentally  found  values  for  the  heat-transfer 
coefficients  were  comparative  with  the  Nusselt  theory. 

Based  upon  research  done  currently,  it  is  evident  that 
present  day  smooth-tube  steam  condensers,  operating  under 
typical  conditions,  have  limitations  in  their  thermal 
efficiency,  due  to  a  large  thermal  resistance  which  occurs 
on  the  tube  side  of  the  condenser.  This  resistance  is 
generally  larger  than  any  of  the  thermal  resistances  that 
occur  on  the  steam  side,  in. the  tube  wall,  those  due  to 
fouling,  or  due  to  noncondensable  gases.  However,  employing 
enhanced  tubes,  the  inside  heat-transfer  coefficient  can  be 
increased  by  100  to  200  percent  over  the  smooth-tube  case. 
The  outside  heat-transfer  coefficient,  on  the  other  hand,  is 
increased  by  only  10  to  50  percent.  In  this  situation,  the 
thermal  resistances  on  the  inside  and  outside  of  the  tube 
can  be  approximately  equal. 


Webb  [Ref.  4]  reported  that  the  dominant  thermal  resis¬ 
tance  in  film  condensation  is  that  of  conduction  across  the 
condensate  film  and,  therefore,  a  surface  geometry  that 
promotes  reduced  film  thickness  will  provide  enhancement. 

Thomas,  et  al,  [Ref.  14]  tested  ammonia  condensation  on 
a  smooth  tube  with  a  wire  wrapped  in  a  helical  manner.  The 
measured  condensing  coefficient  was  approximately  three 
times  that  predicted  by  the  Nusselt  equation  for  a  smooth 
tube.  Surface-tension  forces  draw  the  condensate  to  the 
base  of  the  wires,  which  act  as  condensate  run  off  channels. 

Webb  [Ref.  4]  stated  that,  when  noncondensables  are 
present,  an  additional  thermal  resistance  is  introduced  in 
the  gas  at  the  vapor-liquid  interface.  Mixing  in  the  gas 
film  will  substantially  reduce  this  thermal  resistance. 
Therefore,  the  maintenance  of  high  vapor  velocities,  or 
special  surface  geometries  that  promote  a  higher  heat- 
transfer  coefficient  in  the  gas  film  will  substantially 
alleviate  the  performance  deterioration  due  to 
noncondensables . 

Cunningham  [Ref.  15]  presented  in  his  paper  that,  for 
the  roped  tubes  on  the  vapor  (shell)  side,  the  enhancement 
is  achieved  by  improved  condensate  drainage,  while  on  the 
coolant  (tube)  side  the  helical  ridges  increased  turbulence 
and,  as  a  result,  the  inside  convective  coefficient. 
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Improvements  on  the  condensing  side  up  to  100  percent 
have  been  reported  for  single-tube  tests  [Refs.  16;17]  . 
Although  titanium  has  a  low  thermal  conductivity,  it  pro¬ 
vides  a  high  resistance  to  erosion  and  water-side  fouling. 
Titanium  tubes  with  enhancement  both  inside  and  outside  are 
commercially  available  through  the  Wolverine  Tube  Division 
of  Universal  Oil  Products,  Inc.  The  applications  of  these 
tubes  having  all  the  inherited  properties  of  titanium  are 
promising  for  naval  condensers. 

The  goals  of  this  thesis  were  therefore  to: 

1.  Obtain  baseline  heat-transfer  performance  data  for 
the  test  condenser  utilizing  16  mm  O.D.  smooth 
titanium  tubes. 

2.  Conduct  steam  condensation  tests  with  the  following 
enhanced  tube  geometries  to  determine  steam-side 
heat-transfer  coefficients  in  relation  to  smooth- 
tube  performance: 

a.  Wolverine  "roped"  tubes. 

b.  Wolverine  "roped"  wrapped  with  titanium  wire. 

c.  Smooth  tubes  wrapped  with  titanium  wire. 


The  combined  effect  of  vapor  shear  and  inundation  on  the 
condensate  film  heat-transfer  coefficient  for  cylindrical, 
horizontal  tubes  within  tube  bundles  is  a  very  complex  and 
still  insufficiently-understood  subject,  which  is  of  impor¬ 
tance  to  the  efficient  design  of  steam  condensers.  Although 
many  researchers  have  studied  this  subject  both  theoretic¬ 
ally  and  experimentally,  there  is  no  accurate  methodology 
available  for  predicting  the  condensate-film,  heat-transfer 
coefficient  within  tube  bundles. 

In  1916.  Nusselt  conducted  his  pioneering  analysis  for 
the  simple  case  of  condensation  occurring  on  the  outside  of 
a  single,  isolated,  horizontal  tube.  He  idealized  the 
problem  by  making  the  following  assumptions  for  single  tubes 
as  stated  by  Nobbs  [Ref.  18]. 

1.  The  wall  temperature  is  constant. 

2.  The  flow  is  laminar  in  the  condensate  film. 

3.  Heat  transfer  in  the  condensate  is  by  conduction, 
and  subcooling  may  be  neglected. 

4.  The  fluid  properties  are  constant  within  the 
condensate  film. 

5.  The  forces  due  to  hydrostatic  pressure,  surface 
tension,  inertia,  and  vapor-liquid  interfacial 
shear  are  negligible  when  compared  to  the  viscous 
and  gravitational  forces. 

6.  The  surrounding  steam  and  vapor/liquid  interface 
are  at  saturation  temperature. 


7.  The  film  thickness  is  small  when  compared  with 

normal  tube  diameters  and  the  effects  of  curvature 
are  small. 

Based  on  the  above  assumptions,  Nusselt  predicted  the  famous 
relationship  for  the  heat-transfer  coefficient: 
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In  order  to  simulate  and  analyze  a  tube  bundle, 
Eissenberg  [Ref.  11]  stated  the  following  additional 
assumptions: 

8.  Condensate  drains  as  a  laminar  sheet  from  a  tube  on 
to  the  tube  directly  underneath  in  such  a  way  that 
velocity  and  temperature  gradients  are  not  lost  in 
the  fall  between  tubes. 

9.  The  saturation  temperature  and  the  tube-wall 
temperature  are  constant  for  all  tubes  in  the  bank. 

Jakob  [Ref.  19]  extended  the  Nusselt  analysis  for 
filmwise  condensation  heat  transfer  on  a  vertical  in-line 


row  of  horizontal  tubes  as  shown  in  Figure  la. 

The  above-mentioned  assumptions  were  combined  with  the 


assumption  of  constant  temperature  drop  across  the 
condensate  film  for  all  the  tubes,  and  the  average 
coefficient  for  a  vertical  row  of  N  tubes  was  predicted  to 


Upon,  dividing  equation  (2,2)  by  equation  (2.1),  the  Nusselt 
theory  can  be  expressed  as: 


-1/4  , 
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Equation  (2.3)  can  be  also  expressed  in  terms  of  the 
local  coefficient  for  the  N-th  tube: 


is  -  N  3/4  -  (N-l)  3/4  (2.4) 
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In  reality,  condensate  does  not  drop  off  in  a  continuous 
laminar  sheet,  but  drops  off  instead  by  discrete  droplets  of 
liquid,  as  shown  in  Figure  lb,  depending  upon  the  surface 
tension  of  the  condensate.  These  droplets  create  ripples  in 
the  condensate  film,  and  thereby  decrease  the  performance 
degradation  due  to  inundation.  Based  on  his  research,  Kern 
[Ref.  20]  proposed  a  less  conservative  relationship: 


(2.5) 


or,  in  terms  of  the  local  coefficient  for  the  N-th  tube 


Chen  [Ref.  21]  considered  the  following  conditions: 


1.  the  momentum  gain  of  the  falling  condensate  between 
tubes,  and 

2.  the  condensation  of  vapor  on  the  condensate  between 
tubes, 

and  concluded  that: 
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and 


The  above  approximate  expression,  due  to  Chen,  is  valid  for 
most  ordinary  applications. 

Experimental  work  doen  by  Eissenberg  [Ref.  11] ,  in  order 
to  investigate  the  effects  of  steam  velocity,  condensate 
inundation,  and  noncondensable  gases  on  the  heat-transfer 
coefficient,  revealed  that  condensate  does  not  always  drain 
onto  tubes  aligned  vertically,  but  can  be  diverted  sideways, 
caused  by  local,  vapor-flow  conditions.  The  condensate  thus 
follows  a  staggered  path  as  shown  in  Figure  2. 

Eissenberg,  making  the  assumption  that  the  flow  is 
dominated  by  gravity,  stated  in  this  side-drainage  model, 
that  condensate  strikes  the  lower  tubes  on  their  sides 
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rather  than  their  tops.  Therefore,  the  inundation  effects 
influence  the  condensate  flow  only  on  the  lower  half  of  the 
tubes,  which  transfer  less  heat  than  the  upper  half.  Based 
on  the  above-mentioned  model,  Eissenberg  obtained  the 
following  formula: 

K  _i/4 

r— -  =  0.60  +  0.42  N  (2.8) 

^Nu 

Much  experimental  research  has  been  conducted  for 
studying  the  effect  of  condensate  inundation.  In  general, 
the  obtained  data  are  highly  scattered  as  shown  in  Figure  3. 
Berman  [Ref.  22]  made  a  comprehensive  comparison  of  filmwise 
condensation  data  on  bundles  of  horizontal  tubes,  and  he 
concluded  that  the  wide  variation  in  expermental  data  for 
tube  bundle  inundation  is  caused  by  the  following  variables: 

1.  bundle  geometry  (in-line  or  staggered), 

2.  tube  spacing, 

3.  type  of  condensing  fluid, 

4.  operating  pressure, 

5.  heat  flux,  and 

6.  local  vapor  velocity. 

In  addition,  noncondensable  gases,  and  insufficient  steam 
for  lower  tubes  can  cause  such  scattering  of  data. 


Inundation  Studies 


Nobbs  [Ref.  18]  used  one  active  tube  in  a  dummy  tube 
bundle.  He  simulated  additional  condensate  by  using  three 
porous  tubes.  Based  on  his  results,  he  concluded  the 
following: 

1.  Vapor  velocity  increases  the  condensate  heat- 
transfer  coefficient  on  both  inundated  and 
uninundated  tubes  in  a  given  tube  bundle. 

2.  The  effect  of  inundation  is  to  reduce  the  heat- 
transfer  coefficient. 

3.  The  condensate  drainage  path  is  often  not  vertically 
downwards  but  in  a  diagonal  direction.  This  can 
result  in  tubes  receiving  different  amounts  of 
inundation. 

Marto  and  Nunn  [Ref.  23]  made  a  comprehensive  survey  of 
the  effects  of  noncondensable  gases.  They  stated  that  these 
can  be  classified  into  one  of  two  categories: 

1.  the  introduction  of  an  additional  local  thermal 
resistance  due  to  the  propagation  of  noncondensable 
gases  towards  a  condensing  tube  surface  under  the 
influence  of  a  gradient  in  partial  pressure,  and 

2.  the  cumulative  effect  of  gas  blanketing  where  uneven 
rates  of  condensation  in  a  condenser  bundle  eventu¬ 
ally  lead  to  regions  where  tubes  are  inoperative  in 
a  condensing  role. 

Experiments  done  recently  by  Cunningham  [Ref.  15]  revealed 
that  noncondensable  gases  have  a  less  effect  on  the  perform¬ 
ance  of  finned  tubes  compared  to  smooth  tubes.  Based  on  the 
above  discussion,  it  is  clear  that  the  performance  degrada¬ 
tion  due  to  noncondensables  must  be  taken  into  consideration 
in  realistic  designs  of  tube  bundles. 


Another  very  important  factor  in  the  performance  of  a 


condenser  is  the  effect  of  vapor  velocity.  As  noted 
earlier,  the  vapor  shear  plays  a  beneficial  role.  Berman 
and  Tumanov  [Ref.  24]  conducted  experiments  on  a  single 
horizontal  tube  placed  in  a  bank  of  uncooled  neighboring 
tubes.  For  vapor  in  vertical  downflow,  they  found  a 
relation  between  the  vapor  Reynolds  number,  and  the  heat 
flux  as: 
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Eissenberg  [Ref.  25]  has  stated  that  in  designing  exper¬ 
imental  bundles,  the  combined  effects  of  inundation  and 
vapor  shear  are  very  important.  In  fact,  the  use  of  narrow 
condenser  bundles  to  experimentally  study  vapor  shear  and 
inundation  effects  is  preferred  to  simulate  large  conden¬ 
sers.  However,  small  narrow  tube  bundles  can  create  errors, 
due  to  the  following  causes: 

1.  Wall  flow:  condensate  drainage  may  reach  side 
walls; 

2.  Dummy  tubes:  condensate  inundation  may  disperse; 

3.  Vapor  lanes:  vapor  may  bypass  along  the  side  walls; 
and 

4.  Noncondensable  gases:  they  can  affect  condensation 
even  at  low  concentration  in  the  bulk  stream, 
particularly  if  steam  is  recycled. 


III.  EXPERIMENTAL  FACILITY 


A.  TEST  FACILITY 

The  test  facility  shown  in  Figures  4  and  5,  was  designed 
and  built  by  Morrison  [Ref.  13]  and  modified  by  Noftz 
[Ref.  12]  to  simulate  an  active  tube  column  having  up  to  30 
tubes  in  increments  of  five  tubes  deep  (i.e.,  five,  ten, 
fifteen,  etc.).  Some  elements  of  the  original  test  facility 
were  modified  to  allow  for  more  efficient  operation  of  the 
facility. 

A  detailed  description  of  the  components  used  in  the 
test  facility  is  given  in  Reference  12.  Only  a  short 
description  of  these  components  will  be  found  in  this 
report.  Particular  attention,  however,  will  be  focused  on 
the  experimental  tubes.  Calibration  procedures  for 
components  requiring  calibration  are  outlined  by  Reilly 
[Ref.  7]  . 

B.  STEAM  SYSTEM 

The  steam  system  shown  in  Figure  6  was  modified  from 
Noftz' s  initial  design.  From  the  house  supply  line,  steam 
flows  via  a  19  mm  O.D.  stainless  steel  line  through  a  steam 
supply  valve  (MS-3)  to  a  cast  iron  steam  separator.  The 
steam  continues  through  the  system  past  two  Nupro  bellows 
valves,  which  were  used  in  conjunction  with  the  supply  valve 
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Figure  6.  Schematic  of  Steam  System 


to  regulate  the  steam  supply  pressure.  From  these  valves, 
the  steam  flows  into  the  test  condenser  diffuser. 

The  steam  supply  pressure  was  monitored  by  a  pressure 
gage  just  downstream  of  the  steam  supply  valve;  also  a 
compound  gage  just  after  the  Nupro  valves  was  used  to 
monitor  the  pressure  in  this  line. 

The  operator  had  no  control  over  the  state  point, 
quality  or  noncondensable  gas  content  since  house  steam  was 
used.  However,  especially  at  nights  and  during  weekends, 
the  state  point  of  the  steam  at  the  inlet  to  the  test 
condenser  was  found  to  be  nearly  constant. 

C.  TEST  CONDENSER 

The  test  condenser  shown  in  Figure  7  was  unchanged  from 
Noftz's  initial  design.  Steam  enters  via  the  top,  passes 
through  the  transition  piece,  the  vortex  annihilator,  and 
finally  through  the  diffuser.  The  dimensions  of  the  test 
condenser  were  305  mm  X  305  mm  X  79  mm  and  it  was  made  of 
stainless  steel.  These  dimensions  allowed  for  a  maximum  of 
twenty-seven  16  mm  O.D.  tubes  arranged  in  an  in-line  config¬ 
uration  of  three  columns  of  nine  tubes  each.  For  this 
experiment,  the  in-line  configuration  was  used  with  a  middle 
column  of  five  active  tubes  flanked  on  either  side  by  a 
column  of  five  dummy  tubes.  Just  on  top  of  the  upper  active 
tube,  a  perforated,  distilled  v-  .er  supply  tube  was  posi¬ 
tioned,  and  flanked  on  both  sides  by  dummy  tubes.  In  order 
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to  conduct  the  experiments,  a  square,  in-line  arrangement  of 
tubes,  with  a  pitch-to-diameter  ratio  of  1.5  was  used. 

A  vertical  slot  (Figure  7)  along  the  centerline  of  each 
condenser  end  plate  was  used  for  active  and  perforated  tube 
installation.  The  tubes  were  positioned  using  nylon  tube 
sheets  that  were  attached  to  the  exterior  of  the  condenser 
end  plates. 

To  minimize  heat  losses,  and  also  to  prevent  leaks  from 
the  tube  sheets,  each  test  condenser  side  was  provided  with 
a  nylon  tube  sheet  one-inch  thick  with  six  holes 
(S/D  *  1.5),  having  about  0.5  mm  tube  clearance,  which 
allowed  the  tubes  to  be  easily  slid  into  the  test  condenser 
through  the  tube  sheets.  The  exterior  side  of  each  tube 
sheet  had  grooves  to  support  O-rings.  The  aluminum  tube 
sheets  had  six  holes,  having  a  tube  clearance  of  about 
0.5  mm  and  were  used  as  sealing  plates.  The  diffuser, 
exhaust  plenum,  transition  piece,  vortex  annihilator,  and 
the  exhaust  piping  which  are  shown  in  Figures  8  and  9,  were 
insulated  with  rubber  insulation. 

A  viewing  window  allowed  viewing  of  the  condensation 
process.  A  double-walled  glass  window  was  used,  and  heated 
air  was  fed  through  the  clearance  between  the  two  glasses  to 
eliminate  fogging  on  the  inner  glass. 


Figure  8.  Details  of  Transition  Piece  and  Vortex  Annihilator 
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TO  TEST  CONDENSER  HOTWELL 


Figure  9 


Details  of  Exhaust  and  Condensate  Piping  from 
Exhaust  Plenum 
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D.  TEST  CONDENSER  TUBES 

Four  kinds  of  active  tubes  were  tested  in  this  experi¬ 
ment.  The  tubes  were  manufactured  by  the  Wolverine  Division 
of  Universal  Oil  Products,  and  all  were  made  of  titanium. 

The  first  kind  of  tubes,  shown  In  Figure  10  were  smooth 
titanium  tubes  of  16  mm  O.D.  with  a  1.65  mm  wall  thickness. 
The  second  kind  of  tubes,  shown  in  Figure  11,  were  single¬ 
start,  helically-corrugated  tubes,  designated  as  Low 
Pressure  Drop  (LPD) ,  with  16  mm  O.D.  and  a  1.65  mm  wall 
thickness.  The  third  kind  of  tubes,  shown  in  Figure  12, 
were  also  single-start,  helically-cor rigated  LPD  tubes, 
wrapped  with  titanium  wire  of  1.58  mm  O.D.  The  fourth  kind 
of  tubes,  shown  in  Figure  13,  were  smooth  tubes,  wrapped 
with  1.58  mm  O.D.  titanium  wire  on  the  same  pitch  as  the 
roped  Wolverine  tubes.  In  order  to  wrap  the  titanium  wire, 
the  tubes  were  attached  into  a  lathe,  and  the  wire  was 
welded  at  the  start  of  the  helical  groove.  Then,  under  a 
constant  tension  of  4  kg,  the  wire  was  wrapped  manually  by 
turning  the  chuck  of  the  lathe.  Finally,  the  free  end  of 
the  titanium  wire  was  welded  at  the  end  of  the  helical 
groove. 

The  660-mm-long  active  tubes  were  connected  to  separate 
cooling  water  supply  and  discharge  lines.  The  dummy  tubes 
flanking  the  active  tubes  were  made  of  16  mm  O.D.  stainless 
steel.  These  tubes  served  to  direct  the  steam  flow  so  as  to 
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Figure  11.  Photograph  of  Roped  Titanium  Tube 


y,p~s- 


simulate  actual  conditions  in  a  condenser.  Cooling  water 
was  not  supplied  to  these  tubes  and  they  did  not  penetrate 
the  test  condenser  end  plates.  Special  characteristics  of 
the  Wolverine  tubes  are  listed  in  Table  I. 

E.  PERFORATED  TUBE 

The  perforated  tube  water  supply  system  is  shown  in 
Figure  14.  This  system  consisted  of  a  perforated  tube 
(located  above  the  uppermost  active  tube)  ,  a  water  heater 
which  served  as  a  supply  tank,  a  rotameter  in  order  to 
control  the  amoung  of  water  supplied  to  the  perforated  tube, 
a  pump  driven  by  a  1/2-HP  electric  motor,  a  condensate  pump 
and  associated  piping  system  and  valves.  The  active  length 
of  the  copper-nickel,  perforated  tube  was  305  mm,  which  was 
identical  to  the  length  of  the  test  condenser.  Supply  water 
entered  one  end  of  this  tube;  the  other  end  was  sealed  off. 

The  supply  tank  was  unchanged  from  Noftz's  initial 
design,  except  for  the  addition  of  a  thermocouple  in  order 
to  have  a  direct  indication  of  the  exact  temperature  of  the 
water  supplied  to  the  perforated  tube.  The  water,  heated  to 
the  temperature  of  the  condensate  leaving  the  bottom  tube, 
was  fed  to  the  perforated  tube  by  a  1/2-HP,  electric-motor- 
driven  pump,  via  a  rotameter  and  a  recirculation  valve.  The 
flow  rate  to  the  perforated  tube  was  controlled  by  using  the 
rotometer  and  the  valve  provided  in  the  heater-water 
recirculation  line.  Also,  another  modification  was  made  to 
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DESIGN  CHARACTERISTICS  OF  WOLVERINE 
CORODENSE  TUBES  (TYPE  LPD) 


Outside  diameter  (in) 

Wall  thickness  (in) 

Outside  diameter  (ft) 

Heat  Transfer  Surface  and  Ratio 

a.  Outside  Aq  (ft^/ft) 

*  * 

b.  Outside  to  inside  Aq/A1 
Inside  diameter  D^  (ft) 

2 

Cross  section  for  flow  A  (ft  ) 

cs 

Number  of  groove  starts 

Pitch  (in) 

Groove  radius  (in) 

Depth-transition  length 
(next  to  plain  section)  (in) 


:  5/8 
:  0.035 
:  0.0520 

:  0.163 
:  1.141 
:  0.0435 
:  0.00163 
:  one  +0 
:  0.300 

:  1/32  nominal 


:  1-3/4  max 


Figure  14.  Schematic  of  Perforated  Tube  Water  Supply  System 


feed  the  condensate  collected  in  the  hotwell,  via  the 
condensate  pump,  to  the  supply  tank  (heater)  in  order  to 
facilitate  the  inundation  of  up  to  30  tubes  under  atmos¬ 
pheric  conditions. 

The  main  feature  of  the  perforated  tube  was  to  inundate 
with  condensate  from  above,  in  order  to  simulate  a  tube 
column  of  more  than  five  active  tubes.  When  the  wrapped 
Wolverine  tubes  and  the  smooth  wrapped  tubes  were  tested, 
the  perforated  tube  was  also  wrapped  with  titanium  wire,  by 
using  the  above-mentioned  technique. 

F.  CONDENSATE  SYSTEM 

The  condensate  system  shown  in  Figure  15  was  modified 
from  Noftz's  initial  design,  and  consisted  of  the  test 
condenser  and  hotwell,  the  condensate  pump,  piping  and 
valves. 

The  test  condenser  hotwell  collected  the  condensate 
produced  by  the  test  tubes.  With  valve  C-l  closed,  the 
condensate  mass  flow  rate  from  the  test  condenser  could  be 
measured  using  a  stop  watch.  Opening  the  valve,  the 
condensate  was  fed,  via  the  condensate  pump,  to  the  supply 
tank  for  the  perforated  tube,  or  dumped  into  the  building's 
drainage  system.  The  condensate  line  connecting  the  test 
condenser  and  the  test  condenser  hotwell  were  insulated 
using  rubber  insulation. 


Schematic  of  Condensate  System 


G.  COOLING-WATER  SYSTEM 

The  cooling-water  system  was  a  partially-closed  system 
as  shown  in  Figure  14,  and  it  was  unchanged  from  Noftz's 
initial  design.  House  water  was  used  for  the  test  facility. 
Cooling  water  was  stored  in  a  1.2  meter,  cubical,  plexiglass 
supply  tank,  and  was  pumped  by  a  5-HP,  electrically-driven 
pump  via  51-mm-O.D.  plastic  piping  to  a  manifold.  Five 
rotameters  were  attached  to  the  manifold  to  measure  the  flow 
rate  through  each  active  tube.  Five  regulating  valves  were 
used  to  obtain  any  water  velocity  between  0  and  5  m/s  within 
the  tubes.  The  cooling  water  passing  through  the  rotameters 
was  fed  via  16  mm  O.D.  stainless  steel  tubing,  and  tygon 
tubing  to  the  active  tubes.  The  total  tube  lengths  were 
long  enough  (over  2  meters)  to  ensure  hydraulically 
fully-developed  flow,  and  no  swirling  into  the  tubes.  After 
flowing  through  the  tubes,  the  cooling  water  was  passed 
through  mixing  chambers,  where  the  temperature  profile  was 
destroyed,  to  facilitate  the  steady  measurement  of  the 
outlet  cooling  water  temperatures.  After  the  mixing 
chambers,  the  cooling  water  was  collected  in  the  supply 
tank.  A  7. 5-HP,  electrically-driven  pump  was  used  to  pump 
water  from  the  supply  tank,  through  a  filter,  to  a  cooling 
tower  in  order  to  minimize  temperature  rise  at  the  water 


inlet. 


The  cooling  tower  was  located  outside  the  building  and 
was  composed  of  four  truck  radiators  across  which  air  was 
blown  by  means  of  a  fan.  The  entire  system,  consisting  of 
the  heat  exchanger  and  the  fan  was  enclosed  in  a  wooden 
structure  with  louvered  openings  to  provide  enough 
ventilation. 

The  tygon  tubing  was  secured  to  the  inlet  and  outlet 
sections  of  the  active  tubes  by  means  of  hose  clamps.  The 
outlet  sections,  including  the  mixing  chambers,  were 
insulated  with  rubber  insulation. 

H.  INSTRUMENTATION 

l.  EJLa.w_&afcss 

a.  Foulton  rotameters  were  used  to  measure  the  flow 
rate  of  cooling  water  for  each  active  tube.  Starting  from 
the  top  active  tube  to  the  bottom  one,  the  rotameters  were 
calibrated  giving  100%  maximum  flow  rates  of  66.9+1, 

72.6±1,  71.5+1,  72.8+1,  and  73.3+1  kg/min. 

b.  The  perforated  tube  water  rotameter  was 
calibrated  giving  a  100  percent  volumetric  flow  rate  of 
969±30  ml/min. 

c.  All  rotameters  were  calibrated  using  the 
procedures  noted  in  Appendix  A  of  Reference  7. 


2.  Temperature 


Stainless-steel  sheathed,  copper-constantan 
thermocouples  were  used  as  the  primary  temperature 
monitoring  devices.  Three  thermocouples  were  used  to 
measure  the  inlet  cooling  water  temperature;  and  twelve 
thermocouples  were  utilized  to  measure  the  outlet  cooling 
water  temperature. 

Additionally,  two  thermocouples  were  used  to  measure 
the  steam  saturation  temperature;  one  thermocouple  was  used 
to  monitor  the  condensate  temperature,  and  one  was  used  to 
measure  the  vapor  temperature.  Table  II  lists  the  locations 
monitored. 

A  Gulton  Industries.  West  20,  0-500°  F  temper¬ 
ature  controller  was  used  to  regulate  the  temperature  of  the 
perforated  tube  supply  water.  The  controller  had  a 
manufacturer-stated  accuracy  of  +0.5%  of  the  span  (about 
±1.25°  F) .  It  was  impossible  to  obtain  the  precise 
temperature  for  the  perforated  tube  water,  because  of  the 
low  response  of  the  heating  system. 

During  inundation  runs,  the  condensate  temperature 
rose  with  the  increasing  number  of  tubes.  For  example,  the 
condensate  temperature  after  the  5th  tube  was  94.6°  C 
and  it  increased  to  97.65°  C  after  the  30th  tube. 
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Several  different  types  of  pressure  measurement 
devices  were  used  in  this  facility.  They  were:  a  Bourdon 
tube  pressure  gage  which  was  used  to  measure  the  steam 
supply  pressure,  located  downstream  of  the  steam  supply 
valve;  an  absolute  pressure  transducer  which  was  used  to 
measure  the  test  condenser  pressure;  and,  a  compound 
pressure  gage  was  also  used  to  measure  the  pressure  drop 
downstream  of  the  Nupro  valves.  The  pressure  transducer  was 
calibrated  against  a  mercury  manometer. 

4.  Data  Collection  and  Display 

A  Hewlett-Packard  HP-3054A  Automatic  Data 
Acquisition  System,  with  a  HP-9826  computer  and  a  HP-2671G 
printer  were  used  to  record  and  display  the  thermocouple  and 
pressure  transducer  readings.  The  temperatures  were 
recorded  in  degrees  Celsius,  while  the  pressure  was  recorded 
in  volts  which  was  then  converted,  using  the  calibration 
curve,  into  mm  Hg  absolute.  The  pressure  transducer  was 
assigned  to  channel  019  of  the  HP  Automatic  Data  Acquisition 
System,  while  the  thermocouples  were  assigned  channels  as 


indicated  in  Table  II. 


TABLE  II 


CHANNEL  NUMBERS  FOR  COPPER-CONSTANTAN  THERMOCOUPLES 


Location 


T  . 
Cl 

#1 

T  , 

Cl 

#3 

T  • 
Cl 

#5 

T 

CO 

#1 

T 

CO 

#1 

T  „ 
CO 

#1 

T 

CO 

#2 

T 

CO 

#2 

T 

CO 

#3 

T 

CO 

#3 

T 

CO 

#4 

H3 

O 

O 

#4 

T 

CO 

#4 

T  „ 

CO 

#5 

HJ 

O 

O 

#5 

T 

sat 

T 

sat 

T 

cond 

T 

vap 

Channel 

000 

001 

002 

003 

004 

005 

006 

007 

008 

009 

010 

011 

012 

013 

014 

015 

016 

017 

018 
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IV.  PEQCEDUEES 


A.  INSTALLATION  AND  OPERATING  PROCEDURES 

Prior  to  installation,  the  titanium  tubes  were 
cleaned  using  a  chemical  cleaning  method  [Ref.  26].  The 
steps  in  this  cleaning  process  are  as  follows: 

a.  Swab  the  tube  surface  with  acetone  to  remove  grease. 

b.  Using  a  test  tube  brush,  brush  the  inside  surface  of 
the  tube  with  a  50%  sulfuric  acid  solution  in  order 
to  remove  any  oxides.  Also,  apply  this  solution 

to  the  outside  surface  of  the  tube. 

c.  Rinse  the  inside  and  outside  of  the  tube  with  tap 
water . 

d.  Apply,  using  a  brush,  a  50%  solution  of  sodium 
hydroxide  mixed  with  an  equal  amount  of  ethyl 
alcohol,  at  the  boiling  temperature  {about 
85°  C)  to  the  outside  surface  of  the  tube. 

e.  Rinse  the  tube  with  tap  water. 

f.  Rinse  thoroughly  with  d.i  stilled  water. 

Prior  to  any  run.  the  condenser  tubes  had  to  be  prepared  to 
ensure  filmwise  condensation.  Exterior  and  interior 
surfaces  were  cleaned  to  ensure  proper  wetting  character¬ 
istics.  Also .  the  tubes  were  cleaned  by  running  steam  at 
atmospheric  pressure  through  the  test  condenser  for  about 
twenty  minutes  without  the  cooling  water  running  through  the 
tubes . 
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f  * . 


It  was  also  found  that,  when  drop-wise  condensation 
occurred,  rinsing  the  tubes  using  the  pe rf or ated-tube  water 
was  sufficient  to  restore  film-wise  condensation. 

2.  System  Operation  and  Steadv-State  Conditions 

Complete  operating  instructions  are  listed  in 
Appendix  A.  The  steady-state  condition  was  reached  about 
three  hours  after  initial  system  light-off,  and  about 
fifteen  minutes  after  changes  to  the  cooling  water  or 
perforated  tube  supply  water  flow. 

When  a  steady-state  condition  was  reached,  the  runs 
were  made.  The  duration  of  each  run  was  approximately  one 
minute.  For  each  condition,  for  example,  for  tubes  11 
through  15.  five  consecutive  runs  were  made  and  average 
values  were  computed. 

The  following  data  were  taken  automatically  by  the 
data-acquisition  system: 

a.  the  thermocouple  readings,  and 

b.  the  pressure  transducer  reading. 

Also,  the  following  data  were  read  into  the  computer, 
through  the  keyboard,  for  each  run: 

a.  the  setting  of  each  rotameter,  and 

b.  the  test  condenser  hotwell  levels. 
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DATA  REDUCTION  PROCEDURES 


1. 


by: 


The  heat-transfer  rate  to  the  cooling  water  is  given 


Q  =  mC  (T  -T  . ) 
v  pw  co  ci 


(4.1) 


The  heat-transfer  rate  can  also  be  found  from  the 
overall  heat-transfer  coefficient  by: 


Q  =  U  A  •  LMTD 
O  o 


(4.2) 


where 


LMTD 


(4.3) 


After  combining  equations  (4.1),  (4.2)  and  (4.3),  it 
is  found  that 


(4.4) 


2.  inside..  HeatrTcansfei  Coefficient 

The  heat-tr ansf e r  coefficient  on  the  inside  is 
calculated  from  the  Sieder-Tate  relationship  as  described  in 
Holman  [Ref .  27 ] : 


0.8 


1/3 


0.14 


(4.5) 


N  = 
u 


K.D. 
l  l 


=  Ci 


R 


(JL) 

w 


In  the  above  equation.  is  referred  to  as  the 
Sieder-Tate  coefficient.  The  remainder  of  the  riqht-hand 
side  of  the  above  equation. 


0.14) 


is  referred  to  as  the  Sieder-Tate  parameter,  X. 

For  thermally  and  hydrodynamically  developed  flow  in 
tubes,  the  Sieder-Tate  coefficient  equals  0.027.  However, 
when  short  tubes  are  considered,  as  in  the  case  of  this 
experiment,  fully-developed  conditions  are  not  attained. 
Therefore,  the  value  of  must  be  found  experimentally. 

The  Wilson  plot  was  used  to  arrive  at  the  value  of 
the  Sieder-Tate  coefficient.  The  Wilson  plot  was  developed 
in  1915  [Ref.  281.  It  is  merely  a  plot  of  1/UQ  versus  the 
inverse  of  the  Sieder-Tate  parameter  (which  is  proportional 
to  the  inverse  of  cooling  water  velocity  raised  to  the  0.8 
power) .  The  reasoning  behind  the  Wilson  plot  is  shown  in 
the  following  development. 

Consider  the  equation  for  the  overall  heat-transfer 
coefficient: 
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(4.6) 


For  the  Wilson-plot  method  to  be  successful,  hQ 
must  be  kept  constant.  When  steam  condenses  on  the  shell 
side,  the  above  condition  can  be  achieved  only  if  the  heat 
flux  (q")  is  kept  constant. 

As  required  for  the  Wilson  plot,  when  the  cooling 
water  velocity  increases.  hi  increases.  UQ  increases, 
and  finally  q”  increases - 
Consider 

q"  =  UqlmTD  (4.8) 

In  order  to  keep  q"  constant,  one  must  therefore 
decrease  LMTD  by  lowering  the  steam  saturation  temperature. 
This  requires  tr ial-and-er ror  setting  of  condenser  pressure 
which  is  a  very  difficult  task. 
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To  avoid  this  difficulty  and  yet  arrive  at  a  Sieder 
Tate  constant  not  affected  by  the  varying  q",  a  modified 
Wilson-plot  method  was  used.  This  method  was  developed  by 
Wanniarachchi  [Ref.  291 ,  and  the  required  steps  are  listed 
below: 

1.  Assume  Ci  =  0.027 

2.  Calculate 


LMTD 


-T. 


Cl 


. -T  .  \ 
sat  ci 

-T 

sat  co; 


3.  Calculate 

a  -  *cp  (Tc0  '  Tci> 

4.  Calculate 

U  =  — 2 _ 

uo  Aq • LMTD 

5A.  For  the  first  data  point  only,  do  the  following 

a .  Ass  ume 

c  =  (^)  °-14  =  1 

U 

b.  Calculate  the  Sieder-Tate  parameter 


X 


0.14 


c.  Calculate 


a.  Calculate 


e.  Calculate 


0.14 


f.  Repeat  steps  a  through  e  until  Cf 
assumed  in  step  a  approximately  equals 
cf  calculated  in  step  e. 

g.  Calculate 


(4. 


h.  Assiqn  Q1  =  Q 

The  second  subscript  of  h  on  the  left-hand  side  of 
equation  (4.9)  refers  to  the  first  data  point. 


5B.  For  all  data  points  except  for  first  data 
point,  do  the  following: 

a.  Calculate 


where 


_Q 

Q 


1/3 


b.  Calculate 


c 


Calculate  the  Sieder-Tate  parameter 


X 


R 


e 


0.8  1/3 

r 


0.14 


6.  Repeat  steps  2  through  3  for  all  data  points. 


7. 


Plot  versus 


1 

X 


NOTE:  It  is  more  reasonable  to  plot  h,  versus  X. 

However,  1/h.  versus  1/x  was  plotted  in  order  to 
be  consistent  with  the  original  Wilson-plot  method 


(  -H— '  versus  4  ) 
Uo  X 


8.  Obtain  the  slope  (m) .  of  the  least-squares-f it » 
straight  line. 

9.  Calculate  the  Sieder-Tate  constant. 

D. 

C.  ■  r—^~ 

1  kwm 

10.  Repeat  steps  2  through  9  until  the  assumed  and 
calculated  C^  values  are  approximately  equal. 

A  listing  of  the  computer  program  used  in  this 

method  can  be  found  in  Appendix  D. 

3 .  Outside  Heat-Transfer  Coefficient 

The  heat-transfer  coefficient  on  the  outside  is 

calculated  using  the  following  steps: 

1.  Calculate  the  average  bulk  water  temperature. 

Tb  -  <Tco  +  Tci>/2 
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2.  Evaluate  the  thermophysical  properties  based  on 
the  average  bulk  water  temperature. 

3.  Calculate  the  cooling  water  velocity. 


V 


w 


4.  Calculate  the  water-side  Reynolds  number. 


R 

e 


P.  V  D. 
w  w  1 


w 


5.  Calculate  the  heat  transferred  to  the  cooling 
water. 


Q 


=  m 


T  . ) 
ci 


6.  Calculate  the  heat  flux. 


7.  Calculate  the  Nusselt  coefficient  using  the 
formula: 


h.T 

Nu 


0.651 


1/3 


(4.10) 


NOTE:  The  above  formula  was  employed  since  no  direct 
measurement  of  tube  wall  temperature  was  made. 

To  compare  the  condensate  film  temperature,  as 

required  for  equation  (4.9)  ,  an  iterative  scheme  was  used  as 


outlined  below: 


a.  Assume  T 


C  ■  l  ■!  ■  V  ■  V  9  ^ 


f  ”  Tsat* 

b.  Evaluate  the  relevant  thermophysical  properties, 
which  are  included  in  equation  (4.10)  . 

c.  Calculate  the  Nusselt  coefficient  using 
equation  (4.10)  - 


d.  Evaluate  Tf  using  the  formula: 


'sat 


5 


e.  Repeat  steps  a  through  e  until  Tf  assumed 
in  step  a  approximately  equals  calcu¬ 
lated  in  step  e.  ' 

Calculate  the  inside  heat-transfer  coefficient  using 
the  formula: 


h. 


i 


0.8  1/3 

r 


(4.11) 


In  order  to  determine  the  Cf,  the  average  inner 
wall  temperature  must  be  known.  As  noted  earlier, 
this  temperature  is  not  known  directly  and  it  must  be 
found  iteratively,  as  described  below: 

a.  Assume  a  correction  factor  C^  (say  1)  . 

b.  Calculate  the  h^  using  equation  (4.11) . 

c.  Calculate  the  tube  wall  temperature  using  the 
formula : 


d. 


Calculate  a  new  correction  factor,  at  the 
evaluated  tube  wall  temperature 


/  M , 0.14 


e.  Repeat  steps  a  through  d  until  Cf  assumed 
in  step  a  approximately  equals  calcu¬ 
lated  in  step  d. 


9. 


Calculate  the  log-mean-temperature  difference, 
( LMTD ) . 


LMTD  = 


T  -  T  . 
co  ci 


£n 


T  -T  .  ' 
sat  ci 


(4.12 


T  -T 
sat  c , 


10.  Calculate  the  overall  heat-transfer  coefficient 
using  the  formula: 

Uo  =  LMTD  (4.13 

11.  Calculate  the  outside  heat-transfer  coefficient 
using  the  equation: 


:  The  validity  of  equation  (4.13)  is  based  on  the 
assumption  of  negligible  water-side  fouling 
resistance  and  the  resistance  due  to  noncondensables 

12.  Calculate  the  normalized,  local,  outside  heat- 
transfer  coefficient,  h^h^. 

13.  Calculate  the  normalized  average,  local  outside 
heat-transfer  coefficient  (h.T/h-,). 


C.  DATA- REDUCT ION  PROGRAM 

A  computer  program  was  utilized  to  analyze  the  raw  data 
The  program  was  in  BASIC  language  and  was  run  on  an  HP-9826 
computer  system.  A  peripheral  plotter  was  used  to  plot  the 
results.  The  program  is  presented  in  Appendix  D. 


V.  RESULTS .. AND  DISCUSSION 


A.  SIEDER-TATE  COEFFICIENTS  FOR  SMOOTH  AND  ROPED  TUBES 

The  modified  Wilson  Plots  for  smooth  and  roped  tubes  are 
shown  in  Fiqures  16  and  17.  Of  particular  interest  is  the 
very  small  intercept,  which  represents  the  estimated  fouling 
factor.  For  example,  the  estimated  fouling  factor  for  data 
run  STSD-11  (for  smooth  tubes)  is  1.37X10“®  m^ • K/W 
based  on  the  inside  area.  This  value  is  only  three  percent 
of  the  typical  value  (4.4X10“®  m^.R/W)  used  for 
design  purposes.  This  very  small  value  supports  the  initial 
assumption  of  negligible  fouling  factor  for  the  success  of 
the  modified  Wilson-Plot  method. 

Table  III  is  a  summary  of  the  Sieder-Tate  coefficients 
calculated  for  smooth  and  roped  tubes.  The  average  Sieder- 
Tate  coefficient  was  calculated  to  be  0.029±0.001  for 
smooth  tubes  and  0.061+0.002  for  roped  tubes.  Thus,  the 
roped  tubes  have  a  Sieder-Tate  coefficient  2.1  times  greater 
than  that  for  the  smooth  tubes.  This  increase  is  mainly  due 
to  the  increased  surface  area,  turbulence  and  swirl  effects. 

The  Sieder-Tate  coefficient  derived  for  the  smooth  tubes 
is  7.4  percent  greater  than  the  value  (0.027)  published  in 
the  original,  generalized  correlation.  This  increase  can  be 
easily  explained  by  the  short  condensing  tubes  used  in  this 
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Figure  16.  Modified  Wilson-Plot  for  Smooth  Tubes  (Run  STSD-11) 


TABLE  III 

SUMMARY  OF  SIEDER-TATE  COEFFICIENTS 
FOR  SMOOTH  AND  ROPED  TUBES 


File  Name 


Tube  Tyue 


Sieder-Tate 

Coefficient 


STSD-2 

Smooth 

0.0287 

STSD-10 

Smooth 

0.0291 

STSD-11 

Smooth 

0.0296 

RTSD-3 

Roped 

0.0591 

RTSD-4 

Roped 

0.0589 

RTSD-6 

Roped 

0.0621 

RTSD-7 

Roped 

0.0627 

experiment.  Even  though  the  flow  condition  was  hydro- 
dynamically  fully-developed,  it  was  thermally  developing 
throughout  the  condensing  length,  resulting  in  a  greater 
Sieder-Tate  coefficient. 

It  is  worth  noting  that  the  Sieder-Tate  coefficient 
derived  by  using  the  original.  Wilson-Plot  method  (i.e., 
plotting  1/UQ  versus  1/X) ,  consistently  gave  values  about 
10%  greater  than  the  values  obtained  using  the  modified. 
Wilson-Plot  method. 

B.  SMOOTH  TUBES 

Figure  18  shows  the  variation  of  the  normalized,  local, 
condensing  coefficient  for  five  tubes.  The  data  points  lie 
about  40  percent  above  the  curve  predicted  by  the  Nusselt 
theory.  This  close  agreement  was  considered  to  be  an 
indication  of  the  proper  opepation  of  the  test  apparatus  and 
the  data  reduction  procedures. 

As  discussed  earlier  in  Chapter  II,  the  Nusselt  theory 
for  a  tube  bundle  is  based  on  a  number  of  basic  assumptions. 
For  example,  in  reality,  the  condensate  does  not  fall  as  a 
continuous  sheet  as  assumed  for  the  Nusselt  theory,  but  it 
forms  drops  before  leaving  the  tube.  This  phenomenon  has  a 
tendency  to  result  in  a  larger  condensing  coefficient  than 
that  predicted  by  the  Nusselt  theory. 
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Variation  of  Local  Condensing  Coefficient  with  Tube  Number 
(Run  STNWNI-1) 


Fiqure  19  shows  the  normalized,  average  condensing 
coefficient  for  five  tubes.  These  data  points  show  a  much 
smoother  trend  than  that  for  local  values.  Further,  these 
data  points  are  closer  to  the  Eissenberg  correlation  than 
the  Nusselt  prediction. 

Figure  20  shows  the  variation  of  the  normalized,  local, 
condensing  coefficient  for  30  tubes-  The  data  points  lie  up 
to  35  percent  above  the  curve  predicted  by  the  Nusselt 
theory.  The  outside  heat-transfer  coefficielt  for  the  first 
tube  was  ten  percent  smaller  than  the  value  predicted  by  the 
Nusselt  theory. 

Fiqure  21  displays  the  normalized,  average  condensing 
coefficient  for  30  tubes.  Again,  these  data  points  show  a 
smoother  trend  than  that  for  local  values.  These  points  lie 
between  the  Eissenberg  and  Nusselt  predictions. 

Figure  22  shows  the  least-squares-f it  curve  for  the  data 
points  for  smooth  tubes.  The  exponent  derived  is  -0.154, 
and  it  is  in  close  agreement  with  the  value  of  -0.14  derived 
by  Noftz  [Ref.  12]. 

C.  ROPED  TUBES 

Figure  23  displays  the  variation  of  the  normalized, 
local,  condensing  coefficient  for  five  tubes.  The  data 
points  lie  up  to  15  percent  above  the  curve  representing  the 
Nusselt  theory- 
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Figure  20.  Variation  of  Local  Condensing  Coefficient  with  Tube  Number 
(Run  STNWI-1) 


Figure  21.  Variation  of  Average  Condensing  Coefficient  with  Tube  Number 
(Run  STNWI-1) 


Figure  22.  Least-Squares-Curve  Fit  for  Data  (Run  STNW-1 ) 


Figure  23.  Variation  of  Local  Condensing  Coefficient  with  Tube  Number 
(Run  RTNWNI-3) 


Figure  24  shows  the  normalised,  average,  condensing 
coefficient  for  five  tubes.  These  points  show  a  trend  very 
close  to  the  Eissenberg  relationship. 

Figure  25  represents  the  plot  of  the  normalized,  local, 
condensing  coefficient  under  inundation  conditions  up  to  30 
tubes. 

Figure  26  shows  the  normalized,  average,  condensing 
coefficient  for  30  tubes.  These  data  points  lie  about 
midway  between  the  Eissenberg  and  Nusselt  predictions. 

Figure  27  shows  the  least-squar es-f it  curve  for  the  data 
points  for  roped  tubes.  The  exponent  calculated  is  -0.183. 
Since  the  interest  is  in  the  large  tube  bundles,  the  curve 
fit  was  generated  only  for  tubes  11  through  30.  This  is 
evident  from  the  figure  as  the  curve  fit  shows  poor  agree¬ 
ment  when  compared  to  the  data  points  for  the  first  ten 
tubes. 

D.  ROPED  TUBES  WRAPPED  WITH  WIRE 

The  outside,  heat-transfer  coefficient  for  the  first 
tube  was  11.3  percent  greater  than  the  value  predicted  by 
the  Nusselt  theory.  This  increase  is  in  agreement  with  the 
manufacturer’s  claim  that  the  special  geometry  of  the  roped 
tubes  has  a  thinning  effect  on  the  condensate  film,  result¬ 
ing  in  a  larger  condensing  coefficient. 
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Figure  24.  Variation  of  Average  Condensing  Coefficient  with  Tube  Number 
(Run  RTNWNI-3) 


Variation  of  Local  Condensing  Coefficient  with  Tube  Number 
(Run  RTNWNI-3) 


Figure  26.  Variation  of  Average  Condensing  Coefficient  with  Tube  Number 
(Run  RTNWI-1) 


Figure  27.  Least-Squares-Curve  Fit  for  Data  (Run  RTNWI-1) 


Figure  28  shows  the  variation  of  the  normalized,  local, 
condensing  coefficient  for  five  tubes.  The  data  points  lie 
up  to  63  percent  above  the  curve  predicted  by  the  Nusselt 
theory. 

Figure  29  displays  the  normalized,  average,  condensing 
coefficient.  The  data  points  are  well  above  the  Eissenberg 
correlation. 

Figure  30  shows  the  normalized,  local,  condensing 
coefficient  for  a  bundle  of  30  tubes.  The  data  points  are 
scattered  within  the  limits  of  uncertainty,  as  predicted  by 
the  error  analysis  (see  Appendix  C) . 

Fiqure  31  shows  the  variation  of  the  data  points 
representing  the  normalized,  average,  condensinq  coefficient 
under  inundation  up  to  30  tubes.  The  data  points  lie  up  to 
100  percent  above  the  curve  predicted  by  the  Nusselt  theory. 

Figure  32  represents  the  least-squares-f it  curve,  which 
has  a  derived  exponent  of  -0.039  .  The  c.urve  fit  is  in  good 
agreement  with  all  the  data  points. 

E.  SMOOTH  TUBES  WRAPPED  WITH  WIRE 

The  outside  heat-transfer  coefficient  for  the  first  tube 
was  up  to  six  percent  greater  than  the  value  predicted  by 
the  Nusselt  theory.  This  increase  is  caused  by  the  thinning 
effect  on  the  condensate  film,  resulting  from  the  surface- 
tension  forces  acting  toward  the  wrapped  wire. 


Figure  29.  Variation  of  Average  Condensing  Coefficient  with  Tube  Number 
(Run  RTWNI-3) 


Figure  30.  Variation  of  Local  Condensing  Coefficient  with  Tube  Number 
(Run  RTWI-3) 


Figure  33  shows  the  variation  of  the  normalized,  local, 
condensing  coefficient  for  five  tubes.  The  data  points  lie 
up  to  76  percent  above  the  curve  predicted  by  the  Nusselt 
theory. 

Figure  34  shows  the  normalized,  average  condensing 
coefficient  for  five  tubes.  These  data  points  are  well 
above  the  Eissenberg  correlation- 

Figure  35  shows  the  variation  of  the  normalized,  local, 
condensing  coefficient  for  30  tubes.  The  data  points  lie  up 
to  107  percent  above  the  curve  predicted  by  the  Nusselt 
theory. 

Figure  36  displays  the  normalized,  average,  condensing 
coefficient  for  30  tubes.  The  data  points  lie  above  the 
curve  representing  the  Eissenberg  correlation. 

Figure  37  shows  that  the  curve  fit  is  in  good  agreement 
with  all  the  data  points.  The  derived  exponent  is  -0.037. 

F.  OBSERVATIONS 

l.  Smooth-Tubes 

During  all  runs,  complete  film-wise  condensation  was 
observed  without  any  visible  evidence  of  drop-wise  conden¬ 
sation.  This  was  mainly  achieved  by  the  tube  cleaning 
procedures  used. 

A  slow  rate  of  condensate  droplet  migration,  from 
cooling  water  outlet  end  to  the  inlet  end,  was  observed  at 
the  bottom  of  each  active  tube-  This  problem  was  minimized 


Variation  of  Local  Condensing  Coefficient  with  Tube  Number 
(Run  STWNI-2) 


ariation  of  Average  Condensing  Coefficient  with  Tube  Number 
Run  STWNI-2 ) 


Figure  37.  Least-Squares-Curve  Fit  for  Data  (Run  STWI-1 ) 


TABLE  IV 


SUMMARY  OF  RUNS  WITH  INUNDATION 


File  Name 


Tube  Type 


External  Wire 
Wrapped  or  Not 
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TABLE  V 

SUMMARY  OF 

RUNS  WITHOUT 

INUNDATION 

External  Wire 

File  Name 

Tube  Type 

Wrapped  or  Not 

STNWEV-2 

Smooth 

STNWNI-1 

Smooth 

STNWNI-2 

Smooth 

STNWNI-3 

Smooth 

RTNWNI-3 

Roped 

RTWNI-1 

Roped 

RTWNI-2 

Roped 

RTWNI-3 

Roped 

STWNI-2 

Smooth 

STWNI-4 

Smooth 

STWNI-5 

Smooth 
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TABLE  X 


COMPARISON 

OF  hjj  for  uninundated  tube 

RUNS 

Tribe  Type 

External  Wire 

Wrap  or  not 

K 

Smooth 

No 

0.7986 

Roped 

No 

0.9019 

Roped 

Yes 

0.9651 

Smooth 


Yes 


0.9339 


TABLE  XI 


COMPARISON  OF 


V 


FOR  INUNDATION  TUBE  RUNS 


Tube  Type 


External  Wire 
Wrap  or  not 


Smooth 

Roped 

Roped 


NO 

No 

Yes 


Smooth 


Yes 


TABLE  XII 


EXPONENTS  OF  THE  LEAST-SQUARES-FIT 


File  Name 

Tube  Type 

External  Wire 
Wrap  or  not 

Exponent 

STNWI-1 

Smooth 

No 

0.154 

RTNWI-1 

Roped 

No 

0.183 

-2 

Roped 

No 

0.191 

-3 

Roped 

No 

0.179 

-4 

Roped 

No 

0.185 

RTWI-1 

Roped 

Yes 

0.039 

n 

— 

Roped 

Yes 

0.039 

STWI-1 

Smooth 

Yes 

0.037 

-2 

Smooth 

Yes 

0.056 
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after  leveling  the  tubes,  by  adjusting  the  leveling  nuts  on 
the  test  condenser  support  bracket-  When  the  condensation 
or  the  inundation  rate  increased,  it  was  observed  that  the 
drops  were  formed  at  more  sites  along  the  tubes,  but  the 
droplet  size  for  each  tube  was  nearly  the  same. 

Aqain,  during  the  experiments,  there  was  no  evidence 
of  drop-wise  condensation.  It  was  observed  that  the  falling 
drops  were  coalesced  on  the  lower  surface  of  the  tube, 
especially  in  the  space  between  the  two  successive  grooves, 
and  then  rivulets  were  formed  and  fell  on  the  tube  below. 
Under  condensate  inundation,  the  phenomenon  of  the  rivulet 
formation  was  more  intense.  It  is  also  worth  noting  that 
the  droplet-formation  frequency  was  higher  near  the  cooling 
water  inlet  end.  This  can  be  easily  explained  by  the 
larger,  local,  Sieder-Tate  coefficient  and  by  the  larger 
temperature  difference,  T  -  T  at  the  inlet  compar- 
ed  to  T  -  T  ,  at  the  outlet  end. 

S»  u  CO 

3  -  Roped  Tubes  Wrapped  with  Wire 

The  condensate  was  formed  in  the  space  between  the 
successive  grooves  wrapped  with  wire.  Then  due  to  the 
surface  tension  forces,  the  condensate  was  drawn  to  the  base 
of  the  wire,  and  drops  were  formed  at  the  bottom  surface  of 


the  tubes. 


It  was  observed  that  the  condensate  forms  between 
two  successive,  helically-wrapped  wires,  and  was  then  drawn 
immediately  towards  the  space  between  the  wire  and  the  tube 
surface.  Rivulets  were  drawn  from  the  base  of  the  wire  to 
the  next  tube  below.  It  was  observed  that  there  was  no  drop 
migration  along  the  tube. 

Under  condensate  inundation,  the  above  mentioned 
drop  formation  was  more  intensive,  but  again  no  splashing 
was  observed. 


VI .  CONCLUSIONS 


1.  The  average,  outside,  heat-transfer  coefficient  for 
30  smooth  tubes  was  59  percent  of  the  Nusselt  coefficient 
calculated  for  the  first  tube  in  the  bank. 

2.  The  average,  outside,  heat-transfer  coefficient  for 
30  smooth  tubes  wrapped  with  wire  was  86  percent  of  the 
Nusselt  coefficient  calculated  for  the  first  tube  in  the 
bank. 

3.  The  average,  outside,  heat-transfer  coefficient  for 
30  roped  tubes  was  63  percent  of  the  Nusselt  coefficient 
calculated  for  the  first  tube  in  the  bank. 

4.  The  average,  outside,  heat-transfer  coefficient  for 
30  roped  tubes,  wrapped  with  wire  was  86  percent  of  the 
Nusselt  coefficient  calculated  for  the  first  tube  in  the 
bank . 

5.  The  Sieder-Tate  coefficient  for  roped  tubes  was  2.1 
times  greater  than  that  for  smooth  tubes. 

6.  Wire  wrapping  considerably  improves  the  average 
condensing  coefficient  for  both  smooth  and  roped  tubes  in 
tube  bundles. 

7.  Of  all  cases  investigated  in  this  study,  roped  tubes 
with  wire  wrap  would  be  the  best  candidate  for  designing 
compact  condensers.  However,  the  water  side  pumping  power, 
which  was  not  investigated  in  this  work,  could  be  consider¬ 
ably  higher  than  that  for  the  smooth-tube  case. 
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VII.  RECOMMENDATIONS 


A.  TEST  APPARATUS  MODIFICATIONS 

The  following  test  apparatus  modifications  are 
considered  advisable: 

1.  Redesign  the  test  condenser  flanges  to  allow  a  large 
sealing  area,  and  also  to  allow  for  the  possibility 
of  further  experiments  with  different  S/D  ratios. 

2.  Redesign  the  existing  test  condenser  hotwell  to 
allow  more  reliable  and  convenient  measurement  of 
condensate  flow  rate. 

3.  Install  a  larger  heater  for  the  perforated-tube 
water  supply  tank  to  facilitate  the  simulation  of 
larger  tube  bundles.  Also,  modify  the  temperature 
control  system  for  the  pe rf or ated-tube  water  supply 
tank,  to  allow  for  more  rapid  heating  and  cooling  to 
reduce  long  delays  between  runs. 

4.  Install  a  glass  window  on  the  rear  side  of  the  test 
condenser,  in  order  for  the  operator  to  have  a  view 
on  both  sides  of  the  test  tubes. 

5.  Make  the  system  vacuum  tight  so  that  data  can  be 
taken  at  vacuum  conditions. 


B.  ADDITIONAL  TESTS 

The  following  additional  tests  would  be  important  in 
this  continued  investigation. 

1.  Conduct  tests  with  enhanced  tubes  manufactured  by 
Yorkshire  Imperial  Metals,  Ltd. 

2.  Take  movies  of  the  condensation  process  so  that 
further  conclusions  could  be  drawn  with  regard  to 
condensate  drop  phenomena  and  their  relationship  to 
condenser  performance. 

3.  Conduct  tests  varying  the  pitch  of  the  wrapped  wire 
on  smooth  titanium  tubes  and  determine  the  optimum 
wire  pitch. 
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Withdraw  gas  samples  from  the  test  condenser  and 
analyze  cn  the  gas  chromatograph  to  determine  the 
effect  of  noncondensable  gases  on  the  condenser 
performance . 

Conduct  tests  to  investigate  the  effect  of  vapor 
velocity  on  the  heat-transfer  coefficient. 


APPENDIX  A 


OPERATING  PROCEDURES 


INITIAL  PROCEDURES 

1.  Energize  the  main  circuit  breaker  located  in  Power 
Panel  P-2  on  the  wall  to  the  right  of  the  test 
apparatus. 

2.  Energize  the  circuit  breaker  on  the  left  side  of  the 
old  control  board  by  pressing  the  ON  button. 

3.  Energize  the  following  switches  in  the  control 


panel : 

a. 

#1  - 

Perforated  tube  water  supply  tank  (feed 
pump) . 

b. 

#2  - 

Outlets. 

c. 

#3  - 

Hot  water  heater. 

d. 

#4  - 

Condensate  pump. 

e. 

#6  - 

Cooling  tower. 

f. 

#7  - 

Cooling  water  pump. 

4.  Ensure  all  test  apparatus  valves  are  closed. 

5.  Fill  the  perforated  tube  condensate  supply  tank  with 
distilled  water.  Set  the  temperature  controller  for 
the  perforated  tube  water  supply  tank  at  about 

95°  C,  fully  open  the  recirculation  valve.  P-1, 
and  start  the  pump  to  begin  heating  the  water.  The 
controller  will  have  to  be  reset  to  the  proper 
supply  temperature  once  steady-state  conditions  are 
obtained. 

6.  Fill  the  cooling  water  supply  tank.  This  can  be 
done  by  backfilling  with  valves  CK-1  and  CW-4  open. 

7.  Energize  the  data  acquisition  system. 


B.  OPERATION 


1.  House  Steam 

a.  Open  the  main  supply  valve. 

b.  Open  valve  MS-3  until  the  pressure  gage 
indicates  the  desired  steam  supply  pressure. 

c.  Fully  open  MSD-1  to  drain  any  condensate. 

d.  Open  valves  MS-4  and  MS-5  until  the  desired 
steam  supply  pressure  is  obtained  and  re-adjust 
MS-3  as  necessary. 

2.  Condensate  System.  To  collect  the  condensate  in  the 
test  condenser  hotwell,  operate  the  system  with 
valve  C— 1  closed.  After  a  test  run  is  completed, 
open  valve  C-l  to  drain  the  condensate  by  opening 
valve  C-2  to  the  bilges,  or  operating  the  condensate 
pump  in  order  to  fill  the  perforated  tube  supply 
tank . 

3.  Cooling-Water  System 

a.  Open  valves  CW-1 .  CW-2  and  CW-3 . 

b.  Ensure  valves  CK1-1  and  CW-4  are  closed. 

c.  Energize  the  two  cooling  water  pumps. 

d.  Open  valves  CW-5 ,  CW-6 .  CW-7 .  CW-8  and  CW-9  to 
obtain  the  desired  cooling  rates. 

4.  Perforated  Tube  Water  Supply  System 

a.  Once  steady-state  conditions  are  achieved,  reset 
the  temperature  controller  to  the  proper 
inundation  temperature. 

b.  Adjust  the  rotameter  to  the  required  flow  rate 
for  each  run.  The  supply  tank  recirculation 
valve  may  have  to  be  adjusted  to  achieve  the 
desired  flow  rate,  but  should  never  be  fully 
closed  to  avoid  damaging  the  pump  and  to  ensure 
uniform  water  temperature  in  the  tank. 

c.  Refill  the  supply  tank  as  required,  by  using  the 
existing  piping  filling  system,  by  operating  the 
condensate  pump;  or  filling  the  supply  tank  with 
distilled  water  filling  line. 
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5. 


Miscellaneous 


To  maintain  a  clear  test  condenser  window,  open 
valve  A-2  and  then  energize  and  adjust  the  air 
heater  power  supply.  When  securing,  always  turn  off 
the  power  supply  first  and  allow  the  air  heater  to 
cool  before  securing  valve  A-2. 


SECURING  THE  TEST  APPARATUS 

1.  Secure  the  steam  valves  MS-5.  MS-4.  MS-3  and  the 
main  supply  valve. 

2.  Secure  the  air  compressor. 

3.  Secure  the  perforated-tube  water  supply  system  by 
securing  the  pump,  temperature  controller,  and 
valves  P-1  and  P-4.  Drain  the  supply  tank  by  open¬ 
ing  valve  P-3  (if  desired). 

4.  Secure  the  test  condenser  viewing  window  air  heater 
as  prescribed  above. 

5.  Secure  the  data  acquisition  system. 

6.  Allow  the  test  condenser  to  cool  down  for  about  30 
minutes,  then  secure  the  cooling  water  pumps  and 
close  valves  CS-1 ,  CW-2 ,  CW-3 .  CW-5 .  CW-6 ,  CW-7 , 
CW-8,  CW-9  and  C-4. 

7.  Drain  the  test  condenser  hotwell. 

8.  Secure  all  circuit  breakers. 

9.  Drain  the  cooling  water  system  piping  and  rotameters 
by  opening  valve  CW-10  and  leaving  open  the  cooling 
water  rotameter  supply  valves. 

10.  Drain  the  cooling  water  supply  tank  by  opening  the 
drain  valve  via  the  remote  operating  rod. 

11.  Ensure  all  valves  are  secured. 


APPENDIX  B 


SAMPLE  CALCULATIONS 


A.  RUN  STNWNI-1 

A  sample  calculation  is  performed  in  this  section  to 
illustrate  the  solution  procedure  used  in  the  data  reduction 
program  [Ref  29] . 

The  STNWNI-1  run  was  selected  to  perform  this  analysis: 
INPUT  PARAMETERS 


File 

STNWNI-1 

Pressure  condition 

Atmospheric 

Inundation  condition 

5  tubes 

Month,  date  and  time 

05:11:10:08:50 

Run  number 

1 

Tube  number  1 

2 

3 

4  5 

Inlet  temp  (°C)  28.53 

28.54 

28.55 

28.54  28.53 

Outlet  temp  (°C)  32.55 

32.20 

32.17 

32.05  31.83 

Saturation  temperature 

100.24  (°C) 

Degree  of  superheat 

1.88  (°C) 

Condensate  temperature 

92.41  (oc) 

Static  pressure 

767.33  (mm  Hg) 

The  following  calculations  are  further  limited  only  to 
the  first  tube. 


Determination  of  Average  Bulk  Temperature 


T,(l)  =  T  ,  (1)  +  T  (1)  x0.5 

b  ci  co 


Tb (1)  =  (28.53  +  32.55)  0.5 
Tb(l)  =  30 . 54°C 
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Thermophysical  Properties 


P  =  5.329 
r 


P  =  995.2  kg/nr 
y  =  798  x  10"6  N • s/m2 


C  =  4.178  Kj/kg.K 
pw 


k  =  619  x  10 
w 


W/m*  K 


NOTE:  All  properties  are  calculated  at  the  average  bulk  water 

temperature  from  Table  A. 6  p.  782  [Ref  30] 

3 .  Cooling  Water  Mass  Flow  Rate 
m  =  14.22  Kg/min  =  0.237  Kg/s 

4.  Determination  of  Cooling  Water  Velocity 


Vw  P  A^ 


=  (14.22)  60 _ 

w  (995.2)  (1.56  x  10"4) 

V  =1.53  m/s 
w 


Determination  of  Reynolds  Number 


Re  = 


P  V  D. 
w  w  1 


=  (995.2)  (1.53)  (0.141) 

e  798  x  10"6 


Rg  =  26,904 

Determination  of  Heat  Transfer 
Q  =  A  *  (Tco  -  Tci)  .Cpw 

Q  =  (14.22)  (^)  (32.55  -  28.53) 

Q  =  3,980.5  W 
Determination  of  Heat  Flux 


„„  _  3,980.5 _ 

q  MO. 015875)  (0.305) 

q"  =  261,682.2 

in 


,178) 


Determination  of  Nusselt  Coefficient 


281  x  10 


All  properties  are  calculated  at  the  Saturation 
Temperature.  From  Table  A. 6  p.  782  [Ref  30] 


^u  =  °*651 


_3  ^  2  2 

(682x10  J)  •  (957^8)  (2256. 4x10 J)  (9.81) 


(281xl0“6)  (0.015975)  (261682.2 


hNu  =  11492.6  W/m2K 


Determination  of  T 


f  ,c 


T  =  T 
f,c  sat 


0.5 


T  =  100  24  -  ^61685 . 3 
Af,c  XUU'24  11492.6  U,:> 


T^  =  88.  85°C 
I.  /  c 


Thermophysical  Properties 
k  =  673  x  10-3  W/mK 


1 


0  =  963  Kg/m3 
U  =  326  x  10-6  N • s/m2 
hf  =  2.289.5  x  103  J/kg 

All  properties  are  calculated  at  the  film  temperature 
From  Table  A. 6  p.  782  [Ref  30] 


11.  Determination  of  Nusselt  Coefficient 


=  0.651 


=  0.651 


k  3  P  2h 
Kf  ,pf  fg '9.81 

Wf.D0.q" 


1/3 


hNu  =  10,885.3  W/m‘.K 


-3  3  2  3 

(673  x  10  J) *  (963)  (2289.5  x  10J)  (9.81) 


(326  x  10“6)  (0.015875)  (261685.4) 


1/3 


12.  Determination  of  T 


fjC 


T  =  T 
f,c  xsat 


0.5 


Nu 


T  =  100  24  -  261685 . 4  q  c 
Af,c  iUU*^4  10885.3  U,:5 


T  ,  =  88 . 22°C 

f  ,c 


13.  Determination  of  Logarithmic  Meal  Temperature 
Difference  (LMTD) 


LMTD 


Tco  ~  Tci 
An  /  Tsat~Tci\ 
\Tsat  Tco / 


LMTD 


LMTD 


32.55  ~  28.53 
„ I  ,100.24  -  28.53 
'100.24  -  32.55 


69 . 68°C 


14 .  Determination  of  Overall  Heat  -  Transfer  Coefficient 

u  = 

O  LMTD 

„  _  261685.4 

Uo  “  69.68 


I 
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3755.5 


Uo  = 


Determination  of  Inside  Heat-Transfer  Coefficient 


Assume  Cf  =  1.1 

Ci  =  0.029 


h  =  y..c  o  0-8  p  0.333 
l  Di  Ci*Re  ,Pr  -cf 


hi  =  jrfl'fl  (0*029)  (269040*8)  (5.8290,3333 


hi  =  9303.2  W/m  K 


Determination  of  Inner  Wall  Temperature 


=  Tb  + 


m  (261685.4)  (0.015875) 

iw  ju.54  +  (9303.2)  (0.0141) 

T  =  62 . 2°C 
w 

Determination  of  u  at  the  Average  Wall  Temperature 


w  (62 . 2°C)  =  453  x  10  v  N-s/m 
w 


Determination  of  Correction  Factor 


453  x  10 


B.  RUN  STSD-11 

A  sample  calculation  is  performed  in  this  section  to 
illustrate  the  solution  procedure  used  in  the  modified  Wilson 
plot  program  [Ref  29] . 

The  DP-11  run  was  selected  to  perform  this  analysis 
INPUT  PARAMETERS 


File 

Month,  date  and  time 
Data  point 

Steam  Saturation  Temperature 
Inlet  Temperature 
Outlet  Temperature 
Flowmeter  Reading 


DP-11 

04:01:13:57:20 

#1 

100.01°C 
24 . 42°C 
30 .  75°C 
10% 


1.  Determiantion  of  Average  Bulk  Water  Temperature 


NOTE: 

3 

4 


5 


C  =  4.180  Kj/Kg.K 
pw  J  r 

—  6  2 

v  =  857  x  10  N*s/m 

k  =0.613  W/m-K 
w 

All  properties  are  calculated  at  the  average  bulk  water 
temperature  from  Table  A. 6  p.  782  [Ref.  30] 


Cooling  Water  Mass  Flow  Rate 
Mf  =  6.36  kg/min 

Determination  of  Cooling  Water  Velocity 


M, 


V  = 


W  PA. 


V  =  (6.36)  60 

w 


(997)  (1.56  x  10 


-4 


V  =0.68  m/s 
w 


Determination  of  Reynolds  Number 


p  V  D. 
_  w  w  i 

Re  "  — - 

W 


R  = 
e 


(997)  (0.68)  (0.0141) 


857  x  10 


-6 


R  =  11,154 
e 


.  Determination  of  Heat  Transfer 


Q  =  M.  •  (T  -  T  . )  •  C 
v  f  co  ci  pw 

Q  =  (6.36)  (tJ-)  (30.75  -  24.42)  (4180) 

Go 

Q  =  2,805  W 
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Determination  of  Logarithmic  Mean  Temperature 


erence  (LMTD 


Am  sat  ci  \ 
\Tsat-Tco / 


LMTD  = 


30.75  -  24.42 
/100.01  -  24 . 42  \ 
"UOO.Ol  -  30.75/ 


LMTD  =  72 . 38°C 

Determination  of  Overall  Heat  -  Transfer  Coefficient 


o  tt*D  -L-LMTD 
o 


0  =  _ 280  5_  _ 

uo  ff (6 .015875)  (0.305)  (72.38; 


UQ  =  2,547.715 


m  K 


Determination  of  Sieder  -  Tate  Parameter 


x  =  R  -0.8.p  -0.3333 

X  =  (11 , 154-0 ’ 8 )  (5 . 83_0 - 3333 ) 

X  =  0.0003212 

Determination  of  Inside  Heat  -  Transfer  Coefficient 


based  on  assumed  Sieder  -  Tate  Coefficient 


Assume  Cf  =  1.1  and 
C±  =  0.03 


Kw  _  _  0.8  ~  0,333 


h  =  613  *  i°-~3  (0.03)  (11,1540*8)  (5 .8  3°  *  333 )  (1.1) 
l  0.141 

h.  =  4,463.4 
1  m2  *K 

11.  Determination  of  Average  Wall  Temperature 

Tb  ff.D^ 

2805 

Tw  “  Tb  =  <0. 0141)77.305)  ( 4 , 4  5  3  .TT 

T  -  T,  =  46 . 51°C 
w  b 

T  =  46.51  +  T. 
w  b 

T  =  46.51  +  27.59  =  74°C 
w 

12.  Obtain  uw  at  the  Average  Inner  Wall  Temperature 

u  (T  )  *  375  x  10"6  N- s/m2 
'  w' 

13.  Determination  of  Cfc 

{  \  0,14 

Cfc  "  VWw(TwV 

c  .  (W  X  ItT^  °-14 

fc  \375  x  10"6  / 

Cfc  =  1.12267 

14.  Iterate  for  hf 

Assume  Cf  =  1.12  and 

=  0.03 
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125 


X  =  0.0002861 


3433.2 


Uo  = 


Uo  =  3,088.4  w/m  K 


Determination  of  Sieder-Tate  Parameter 


X  =  Re”°*8-Pr~°‘3333 


X  =  (17 , 082-0 ‘ 8)  (5.823-0*3333) 

X  =  0.000228 


Determination  of 


,  s  V3 

i  ±.±(?°)  _B  f° i) 
hi  '  Do  ho  W  w  VV 

57  *  [Wstt  -  mrr  (sttott)  -o.3oooo«2»5  ] 


015875 

0141 


pj-  =  0.0002095 


Determination  of  Average  Wall  Temperature 


T  =  T  +  - — - 

w  ff  •D^*L*h^ 


_  „„  .  3433.2 

Tw  “  ir«  (0-0158  75)  (0.305)  (4,77  3.26] 


T  *  7  4 . 20°C 
w 


Determination  of  at  the  Average  Wall  Temperature 


u  (T  )  =  375.2  x  10 
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B.  UNCERTAINTY  IN  THE  REYNOLD'S  NUMBER 


P  V  D. 
v/w  l 

Re  [x 
Kw 

The  probable  uncertainty  is  given  by: 


SR, 


R. 


SV. 


<£*>2  +  «  V. 


V  ♦  ( 


w 


SD .  -  *  „ 

-D1’  +  '^>2 


1/2 


The  following  uncertainties  were  assigned  to  the  variables 

3 


Sp  =  +3  kg/m' 


5Di  =  ±0.0001  m 
6m  =  ±8xl0-6  N* s/m2 
SVw  =  from  equation  (C.2) 


C.  UNCERTAINTY  IN  HEAT  TRANSFER 


Q  =  m(T  -T  .  )C 

CO  Cl  pw 


The  probable  uncertainty  is  given  by 

5T 


ifi  = 
Q 


co 


v  m;  '  T  -T  . 


SC 


1/2 


)2  *  +  (-T^)2 


T  -T 

CO  Cl  CO  Cl  pw 

The  following  uncertainties  were  assigned  to  the  variables 
5m  =  +0.01  kg/s 


ST  =  +0.025  °C 

CO  — 

ST  .  =  +0.025  °C 
ci  — 

6Cnw  *  ±8  J/Kg*c 


D.  UNCERTAINTY  OF  THE  HEAT-FLUX 


q"  = 


ffD  L 
o 


The  probable  uncertainty  is  given  by: 


1/3 


fig"  _  l 

q"  it 


5D. 


<^)2  +  <^2)2  +  (ii)2 


The  following  uncertainties  were  assigned  to  the  variables 


6D  =  +0.0001  m 
o  — 


6L  =  +0.0001  m 

5Q  =  as  found  from  equation  (C.4) 


E.  UNCERTAINTY  OF  h 


Nu 


h.,  =  0.651 

Nu 


kfp2hfg'9 


1/3 


^fDoq 

The  probable  uncertainty  is  given  by: 


Nu 


=  0.651 


Nu 


(-t1)2  +  (i^)2  +  (  !ufg)2  +  (-4?)2  +  (• 


<SD 


+  (_ I2)2  +  { _A9L!1 )  2 

1  3D  ;  1  3q"; 

o 


3  p 
,  1/2 


3h 


fg 


3g‘ 


The  following  uncertainties  were  assigned  to  the  variables 


5kf  =  +0.0012  w/m*k 
.2 


6  p  =  +3  kg/m 


<5  h 


fg  =  ±0 . 48  J/kg 


5  Hf  =  +8xl0-6  N-S/m2 


+0.0001  m 


<Sg  =  +0.001  m/s 

6q"  =  as  found  from  equation  (C.5) 


F.  UNCERTAINTY  OF  T.  .  , 

1 1 ini 


T,  =  T 
fc  sat 


q" 


Nu 


0.5 


The  probable  uncertainty  is  given  by: 


r  6T  .  n  » 

i-^)2  +  i1^)2  ♦ 

sat  4 


NU)  2 
Nu 


1/2 


The  following  uncertainties  were  assigned  to  the  variables 

STsat  =  ±0'025  °C 

6q"=  as  found  from  equation  (C.5) 

5hNu  ~  as  found  from  equation  (C.6) 


G.  UNCERTAINTY  OF  OVERALL  HEAT-TRANSFER  COEFFICIENT 
q" 

rj  =  —a - 

o  LMTD 

The  probable  uncertainty  is  given  by: 


1/2 

aUQ  _  f/ia’\24.  ,  6  (LMTD)  '  2 

U  '  q";  1  LMTD  ’ 

o  L  J 

The  following  uncertainties  were  assigned  to  the  variables 


5q  =  as  found  from  equation  (C.5) 
6 (LMTD)  =  as  found  from  equation  (C.9) 


H.  UNCERTAINTY  FOR  LOGARITHMIC  MEAN  TEMPERATURE  DIFFERENCE 
( LMTD ) 


LMTD  = 


T  -T 

co  Cl 

T  -T  . 

ln(Ts_^) 

S  CO 


The  probable  uncertainty  is  given  by: 


5  LMTD 
LMTD 


$T  (T  .-T  ) 

S  Cl  CO 


(T  -T  . ) (T  -T  )  T  -T  . 

v  s  ci  s  co  .  ,  s  ci . , 

111  \  m  m  / 


1  T  -T 
S  CO 


ST 


Cl 


T  -T  .  T  -T  ■ 

s  ci  ,  ,  s  ci 

ln( 


T  -T 
s  co 


St 


1/2 


CO 


(T  -T  )  T  -T  . 

>  S  CO  ,  S  Cl,/ 

ln»  _ m  ' 


‘T  -T 
S  CO 


( C .  9  ) 


The  following  uncertainties  were  assigned  to  the  variables: 

6TS  =  +0.025  °C 

ST  .  =  +0.025  ’C 
ci  — 

ST  =  +0.025  °C 
co  — 


UNCERTAINTY  OF  INSIDE  HEAT-TRANSFER  COEFFICIENT 


h. 


..  kw  •  r  p°*8p  0-333,  u  ,  .*• .  1 4 
i  Di  uiKe  r  mw; 
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The  probable  uncertainty  is  given  by: 


.  f<^»2  ♦  <^>  ♦  <^>2  -  <^!£>2  - - 

i  L  Kw  i  e  r  1 

1/2 

0 . 146 ( p/p  )  2  ‘ 

+  <  M/M  — > 

pm  -1 

The  following  uncertainties  were  applied  to  the  variables 


(C.10) 


5k  =  +0.0012  w/m-k 
6D^  =  +0.0001  m 
SC.  =  +0.0001 

5R  =  as  found  from  equation  (C.3) 
e 

6Pr  =  +0.17 

6(p/nw)  =  8xl0~6  Ns/m2 


UNCERTAINTY  IN  TEMPERATURE  DIFFERENCE 


q"  D 
M  o 


6AT 

AT 


^>2  *  <-KT»2  +  < 


;D„  O  5 

TT1  +  ‘“iT1 

o  1 


;c.n ) 


The  following  uncertainties  were  assigned  to  the  variables: 

5q" =  as  found  from  equation  (C.5) 

5h^  =  as  found  from  equation  (C.10) 

5  D  =  +0.0001  m 
o  — 

<5°i  =  +0.0001  m 
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K.  UNCERTAINTY  IN  OUTSIDE  HEAT-TRANSFER  COEFFICIENT 


The  following  uncertainties  were  assigned  to  the  variables: 

SUQ  =  as  found  from  equation  (C.8) 

6h^  =  as  found  from  equation  (C.10) 

<5R  =  +0.00001  m2-k/w 

w  — 

L.  UNCERTAINTIES  FOR  THE  NORMALIZED  LOCAL  HEAT-TRANSFER 

COEFFICIENT  h.,/h . 

N  i 

This  ratio  is  simply  the  heat  transfer  coefficient  of  a 
given  tube,  N,  divided  by  that  of  the  first  tube,  i.e.,  for 
the  fifth  tube,  N=5  and: 
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An  application  of  equation  (C.l)  results  in  the  following 
equation: 


5  <Vhi> 

(W 


(^>2  ♦  ^ 
hi  Ni 


(C.13) 


Note:  Equation  (C.13)  is  valid  only  for  N7/2. 
f  /V.  \  P  fV.  Jji  T  1/2 


For  example: 


6(h,/h.) 


6h0  ^ 

<-e->  +  <Tg» 


M.  UNCERTAINTY  IN  THE  NORMALIZED  AVERAGE  HEAT-TRANSFER 
COEFFICIENT,  h*./!^ 

The  normalized  average  heat-transfer  coefficient  is  obtained 
for  the  Nth  tube  by  taking  the  average  of  the  heat-transfer 
coefficients  of  the  first  N  tubes  and  dividing  this  by  the 
heat-transfer  coefficient  of  the  first  tube: 

(h,  +  h2  +. . .+  hN)N 


Applying  equation  (C.l)  to  the  above,  the  following  equation 
results : 


<Vhi»  Ty  .«hi 

<Vhi>  "  I  i-i  T  hj 

i=l  1 


5hl  2 


where  N  =  the  tube  number.  For  example: 


(h_/h,)  r  5h.  ? 

— — —  =  (*r4c-)  + 


5  h  «  ^  6h.  . 

(r-^4-)2  +  (-TT^r 


For  Run  STNWNI-1 


VW  *  !*54  ±  °.ll  m/s 

Re  =  27546  +  2121 

Q  =  3934  +  169  W 

q"  *  267830  +  11516  w/  2 
“  m 

h*.  =  10802  +  164  W/  2, 

LMTD  =  69.7  +  0.6  °C 

U  =  3843.7  +  165  W/  2, 

o  —  'm*k 

h.  =  9381  +  122  W/  2. 

i  —  'm-k 
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APPENDIX  D 


COMPUTER  PROGRAMS 


DATA  REDUCTION  PROGRAM 


1000*  FILE  NAME :  DRP 

1010!  REVISED:  May  20.  1983 

1020  COM  /Ci/  C < 7 > 

1030  DIM  7ci(2'.Tco<4,2).Ti<4).Mft(4>.Vw<4>.Ho<4) 

1  040  DIM  To<  4 ) .  Ts< 1 > . Tb<  4  > . R3<  4  > .  R4<  4 ) ,S3<  4 > . S4  <  4  > 

1050! 

1080!  ASSIGN  COEFFICIENTS  FOR  THE  3-TH  ORDER 
1070!  POLYNOMIAL  FOR  TYPE-T  (COPPER- CONST ANT AN) 

1080!  THERMOCOUPLES 

1090  DATA  0. 10086091 .25727.94369.-767345.3295.78025595.81 
1100  DATA  -9247486589 .6. 97688E+1 1 , -2.661 92E+1 3 .3 .94078E+1 4 
1110  READ  C< * > 

1120! 

1130!  ASblbN  Full-SCALE  FLOW  RATES  THROUGH  THE  5 
1140!  FLOW  METERS  <ka/run> 

1150  DATA  66.36.73.35.72.44.72.52.72.24 
1160  READ  Mf t(*> 

1170! 

1180!  ASSIGN  SIEDER-TATE  COEFFICIENT  AND  EXPONENT- 
1190!  FOR  REYNOLDS  NUMBER 
1200  Ci * . 023 
1210  Ex-. 8 
1220! 

1230 !  ASSIGN  GEOMETRIC  VARIABLES 
1240  D i - .0*41  !  Inner  diameter  <m> 

1250  Do*. 015875  !  Outer  diameter  <m) 

1260  Ktn=21 .9  !  Thermal  conductivity  of  titanium  <W/m-K> 

1270  L-.3Q5  !  Condensing  length  <m> 

1280  Nc*3.3333  !  Number  of  unit  ceils  across  condenser  width 

1290  Pt-1.5  !  Transverse  tube  p i tch- to-diameter  ratio 

1300! 

1310!  COMPUTE  THE  MINIMUM  S’EAM  FLOW  AREA  IN  THE  TEST  CONDENSER  (m'2) 

1320  AmT »Nc»Do*<Pt-PI/<4»Pt > >*L 

1330! 

1340!  COMPUTE  INSIDE  AREA  AND  WALL  RESISTANCE 

1350  A  i  -PI*D i ' 2/4 

1360  Rw*Do*LGG(Do/Di ) /<2»Ktm) 

1370! 

1380  PRINTER  IS  701 
1390  CLEAR  709 
1400  BEEP 

1410  INPUT  "ENTER  MONTH.  DATE.  AND  TIME  <MM: DD :HH :MM :SS>” . T imeS 
1420  OUTPUT  709 :"TD": Time'S 
(430  BEE? 

1440  INPUT  "ENTER  THE  INPUT  MODE  ( 1 -3054A .2-FILE)” . Im 
1450  IF  Im-2  THEN 
1460  BEEP 

1470  INPUT  "ENTER  THE  NAME  OF  THE  EXISTING  DP'.A  FILE”  .OlddataS 

1480  PRINT  USING  "10X.""This  analysis  was  performed  for  data  stored  in 

10A" :01ddataS 

1490  ASSIGN  5Fiie2  TO  OlddataS 
1500  END  IF 
1510  IF  Im-1  THEN 
1520  BEEP 

1530  INPUT  "GIVE  A  NAME  FOR  THE  DATA  FILE  TO  BE  CREATED" .NewdataS 
1540  CREATE  BOAT  NewdataS. 20 
1550  ASSIGN  jFilel  TO  NewaataS 


file 


1560  END  IF 
1570  BEEP 

1530  INPUT  "GIVE  A  NAME  FOR  THE  OUTPUT  FILE’VFi le_outS 
1590  BEEP 

1600  INPUT  "ENTER  THE  PRESSURE  CONDITION  ( 1 -ATM . 2* VACUUM > “ . Mo 
1610  IF  Mo-1  THEN  PRINT  "  Pressure  condition:  ATMOSPHERIC" 

1620  IF  Mp-0  THEN  PRINT  "  Pressure  condition:  VACUUM" 

1630  BEEP 

1640  INPUT  "ENTER  THE  INUNDATION  CONDITION  <1-5  TUBES.  2-30  TUBES)". Mi 
1650  IF  Mi-2  THEN  PRINT  "  Inundation  condition:  30  TUBES" 

1660  IF  Mi-1  THEN  PRINT  "  Inundation  condition:  5  TUBES" 

1670  CREATE  BDAT  Fiie_outS.6 

1680  ASSIGN  ?File3  TO  Fiie_outS 

1690  Ja-0 

1 700  Nrun-0 

1710  FOR  I -0  TO  4 

1720  S3< I >-0 . 

1 730  $4<Ii-u. 

1740  NEXT  I 
1750  Repeat:  1 
1760  Nrun-Nrun+I 
1770  OUTPUT  709 : "TD" 

1780  ENTER  709:TimeS 
1790  PRINT  "  " 

1800  PRINT  USING  " 1 0X . “"Mon th .  date,  and  tine:  15A" : Time$ 

1810  IF  In-2  THEN  Rdf 
1820  BEEP 

1830  INPUT  "ENTER  FLOW  METER  READINGS  <AS  PERCENTAGES)" .Fnl .Fn2 
1840  IF  Nrun  MOD  5-1  AND  Mi-2  AND  Nrun>5  THEN 
1850  BEEP 

1860  INPUT  "ENTER  FLOW  RATE  FOR  POROUS  TUBE  (AS  A  PERCENT )".Fpt 

1870  OUTPUT  *FiieT:Fpt 

1830  Mpt--9. 361613+10. 076742*Fpt 

1890  END  IF 

1900  DISP  "START  COLLECTING  CONDENSATE" 

1910  BEEP 
1920  WAIT  20 

1930  OUTPUT  709 : "AR  AFO  ALI9" 

1940  OUTPUT  722:"F1  R1  T1  Z1  FL1" 

1950! 

1960!  READ  INLET  WATER  TEMPERATURES 
1970! 

1980  FOR  1=0  TO  2 
1990  OUTPUT  709: "AS  SA" 

2000  ENTER  722:Tci(I) 

2010  CALL  Tvsv<Tci<I>) 

2020  Tci < I ) =FNTemp< Tci ( I ) , I ) 

2030  NEXT  I 
2040! 

2050!  READ  OUTLET  WATER  TEMPERATURES 
2060! 

2070  1 1*2 

2080  FOR  1-0  TO  4 

2090  IF  1-0  OR  1-3  THEN 

2100  Iu-2 

2110  ELSE 

2120  Iu-1 

2130  END  IF 

2140  FOR  J-0  TO  Iu 

2150  Ii-Ii+I 

2160  OUTPUT  709 : "AS  SA" 
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2170  ENTER  722:Tco<I.J> 

2180  CALL  Tvsv<  Tco< I . J> ) 

2190  Tco< I . J>  *FNTemp( Tco( I . J) , I i ) 

2200  NEXT  J 
2210  NEXT  I 
2220! 

2230!  READ  STEAM  TEMPERATURES 
2240! 

2250  FOR  1-15  TO  16 
2260  OUTPUT  709: "AS  SA" 

2270  ENTER  722:Ts<I-15> 

2280  CALL  TvSv< TS(I-15> > 

2290  Ts< I- 1 5) -FNTemp( Ts< I- 15) . II 

2300  NEXT  I 

2310! 

2320!  READ  CONDENSATE  TEMPERATURE 
2330! 

2340  OUTPUT  70S : “AS  SA" 

2350  ENTER  722:Tcon 
2360  CALL  Tvsv(Tcon) 

2370  Tcon-FNTempdcon .  1  7 ) 

2330! 

2390!  READ  VAPOR  TEMPERATURE 
2400! 

2410  OUTPUT  709: "AS  SA" 

2420  ENTER  722: Tv 
2430  CALL  Tvsvdv) 

2440  Tv-FNTenp< Tv, 13) 

2450! 

2460!  READ  VAPOR  PRESSURE 
2470! 

2480  OUTPUT  709: "AS  SA" 

2490  ENTER  722:P_voits 
2500! 

2510!  COMPUTE  AVERAGE  WATER  TEMPERATURES  AT  INLET 
2520! 

2530  Ti<0)-Tci(0) 

2540  Ti(1>-(Tci<0>+Tcl<1>>*.5 

2550  Ti(2)-Tci<l> 

2560  T i < 3) * < Tc i < 1 ) +Tc  i  ( 2> > * .  5 
2570  Ti<4)-Tci<2) 

2580! 

2590!  COMPUTE  AVERAGE  WATER  TEMPERATURES  AT  OUTLET 
2600! 

2610  FOR  1-0  TO  4 
2620  IF  1-0  OR  1-3  THEN 

2630  Tod )  -  <  Tco(1 ,0)+Tco<  1 , 1  >+Tco<  1 .2)  >**  .3333 
2640  ELSE 

2650  To<I)”<Tcod,0)+Tcod,1  )>«.5 
2660  END  IF 
2670  NEXT  I 

2680  Tsa*(Ts(0)+Tsd))«.5 

2690  Pvap-FNPvsv< P_vo 1 ts) 

2700  Tsat-FNTvsp< Pvap ) 

2710  Dsup-Tv-Tsat 
2720! 

2730!  READ  INFORMATION  FOR  CONDENSATE  FLOW  RATE 
2740! 

2750  BEEP 

2760  INPUT  "ENTER  INITIAL  AND  FINAL  LEVELS  IN  HOT  WELL  1".H1.H2 
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2770  Dh-H2-H1 

2780  IF  Nrun  MOD  5-1  THEN  Msum-0 
2790  Mf1-540.4836*Dh 

2800  Md I -Mf ! *FNRhow< Tsat- 10>*1 .0E-6/60 
2810  Msum-Msum+Mf 1 

2820  IF  Mi >2  AND  Nrun<>30  AND  Nrun  MOD  5-0  THEN 
2830  Mave-Msum/5 

2840  Se  t  -  <  Mave«FNRhoui(  T  sat  -  I0)/10'S*-.  03238 )  / .  0421 32 
2850  END  IF 
2860* 

2870  Rdf:  ! 

2880* 

2890  PRINT  USING  "1  OX . ""Run  number  -  "" , DD" :Nrun 

2900  PRINT  "  Tube  *  :  1  2  3  4  5" 

2910  IF  Im-2  THEN 

2920  IF  Nrun  MOD  5-1  AND  Mi-2  AND  Nrun>5  THEN  ENTER  SFile2:Fot 
293U  ENTtR  «*F i ie2 : Ti <  *  > . To(-) . Tsa. Toon , Tv .Pvsp , Tsat ,Dsup .Fn1 .  Fm2 
2940  ENTER  >FiIe2:HI.H2 
2950  END  IF 

2960  IF  Im* 1  THEN  OUTPUT  3F i lei : T i <» ) . To(*> . Tsa. Tcon , Tv .Pvap . Tsat .Dsup .Fml ,Fm2 
2970  PRINT  USING  "1  OX . ""Inlet  temp  (Deg  C>  : ”" . 5< DDD . DD  . 2X ) ”  :  T i < - ) 

2980  PRINT  USING  "1  OX . ""Out  let  temp  (Deg  C> : . 5( DDD . DD . 2X ) " : Tot *> 

2990  PRINT  USING  "1  OX . '"'Saturat ion  temperature  -  "“.3D.DD.""  (Deg  C) . ';Tsat 

3000  PRINT  USING  " 1  OX . ""Degree  of  superheat  -  "".3D.DD.""  (Deg  C> . :Dsup 

3010  PRINT  USING  "1  OX . '"'Condensate  temperature  •  ”",3D.DD,"”  (Deg  C) . :Tcon 

3020  PRINT  USING  "1  OX .""Static  pressure  -  "".3D.DD.""  (mm  Hg> . ;Pvap 

3030  IF  Im-1  THEN  OUTPUT  *Fi lei : HI ,H2 
3040! 

3050*  CALCULATE  AVERAGE  BULK  TEMPERATURES 
3060* 

3070  FOR  1-0  TO  4 

3080  Tb(I>-(Ti<I)+To(I)>».5 

3090  NEXT  I 

3100! 

3110  IF  Mi-1  OR  (Mi-2  AND  Nrun<6>  THEN  As!-0. 

3120  IF  Mi-2  AND  Nrun>5  THEN  Sl-Asl 

3130  Sl-Asl 

3140  FOR  J-0  TO  4 

3150  IF  J-0  THEM  Cuf-Fml 

3160  IF  J-1  THEN  Cwf*Fir.2 

3170  Mf-Mf t( J)-Cuf/( 100-60) 

3180  Tx-Tb(J) 

3190  Vy(J)-Mf/(FNRhou(Tx)*Ai) 

3200! 

3210!  CALCULATE  INSIDE  AND  OUTSIDE  COEFFICIENTS 
3220* 

3230  Rew*FNRhou(  Tx )  **Vu<  J  )»D  i  /FNMuw<  Tx ) 

3240  Cf-1. 

3250  Q-Mf'*FNCpij<Tx)»(To(J)-Ti(J)> 

3260  0p-0/(PI«Do*D 

3270  IF  (Mi-1  OR  (Mi-2  AND  Nrun<6>>  AND  J-0  THEN 
3280  Tfiln-Tsat 
3290  Kf-FNKudf  ilm> 

3300  Rhof -FNRhowt Tf i lm) 

3310  Hf g-FNHf  g<  Tsat ) *  1 000 
3320  Muf -FNMuut  T film) 

3330  Hnu- . 651 «(Kf  '3-Rhot '2«Hfg*9 .81 /(Muf-Do-Qp) )  .3333 

3340  Tf i imc»Tsat-Qp/Hnu- .5 

3350  IF  ABSftTf ilmc-Tf iIm)/Tf ilmc)>.01  THEN 

3360  T  f i lm-Tf i lmc 
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<W/nT2> 


3370  GOTO  32S0 
3330  END  IF 

3390  PRINT  USING  "1  OX .  ""Nussel  t  coefficient  for  first  tube  -  ""  .5D.  L) .  *,M 

. :Hnu 

3400  END  IF 
3410  IF  J-0  THEN 

3420  IF  Hi  *  1  AND  Nrun-1  THEN  Hol*0. 

3430  PRINT  ”  Tube  Uu  Heat  flux  Cond  coef  Rt  R2 

RR‘‘ 

3440  PRINT  "  *  <m/S>  <W/r>'2>  (W/m'2-K)" 

3450  END  IF 
3460  Muw-F'iMuw  <  Tx  ) 

3470  Hi-FNKw<  Tx  >/Di*Ci«Rew  Ex»<FNPru(Tx) ) ' . 3333»Cf 

3480  Dt-Qp/H i*Do/D  i 

3490  Cfc-<Muw/<FNMuu<Tx+Dt>>>'.14 

3500  IF  ABS<  <Cf-Cfc)/Cfc)>.01  THEN 

3510  Cf-<Cf+Cfc>*.5 

352u  guiu  J4/ u 

3530  END  IF 

3540  Lmtd-<  To< J)-Ti<J) >/L0G< < Tsat-Ti ( J) >/<Tsat-To< J) > > 

3550  Uo*Qp/Lnta 

3560  Ho<J)-t ./< 1 ./Uo-Do/<Di*Hi)-Rw> 

3570  Rr-Uo/Ho<J> 

3580  S1-S1-Ho<J> 

3590  IF  Nrun  MOD  5-1  THEN 

3600  IF  Mi - 1  OR  < Mi *2  AND  Nrun-31)  THEN  Ja-0 
3610  IF  Mi*2  AND  5<Nrun  AND  Nrun<30  THEN  Ja-Nrun-1 
3620  IF  Mi-2  AND  35<Nrun  THEN  Ja-Nrun-1 
3630  END  IF 

364C  IF  Mi  *  1  OR  <30<Nrun  AND  Nrun<36  AND  Mi-1)  OR  Nrun<6  THEN 
3650  R1-Ho(J)/Ho<0> 

3660  R2-S! /f  < J+1 +Ja)*Ho<0) ) 

3670  ELSE 

3680  R 1 -Ho<  J ) /Hoi 

3690  R2-S1 /< < J  +  1 +Ja)«Ho1 > 

3700  END  IF 
3710! 

3720!  PRINT  RESULTS 
3730! 

3740  PRINT  USING  "1 1 X . DO , 4X , DD . DD . 2X ,2(D . 5DE .2X> ,3<Z.4D.2X)”; J+1+Ja.Uio<J) .Qp.Ho 

<J) . R 1 .R2.Rr 

3750! 

3760! 

3770! 

3780  IF  Mi -2  AND  Nrun<6  AND  J-0  THEN 

3790  Hoi -Hoi +H'o<  0 ) /5 

3800  END  IF 

3810  FOR  K-0  TO  4 

3820  IF  K-J  ThEN  S3(K)-S3<K)+R1 

3830  IF  K-J  THEN  S4(K)-S4<tO+R2 

3840  NEXT  K 

3850  NEXT  J 

3860  IF  Nrun  MOD  5-0  THEN 
3870  FOR  K-0  TO  4 
3880  R3<K)-S3<K)/5 

3890  R4<K)-S4(K>/5 

3900  S3<K)-0. 

3910  S4<  K ) -0 . 

3920  NEXT  K 

3930  IF  Mi-2  AND  Nrun  MOD  5-0  AND  Nrun<>30  THEN  As! -Nrun»R4<4>-Ho! 

3940! 
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03t>u:  rKiNl  HVtrtHljt  KHfiUb 
3960! 

3970  PRINT  "  " 

3980  PRINT  "  Tube  *  R3  R4" 

3990  FOR  J-i  TO  5 

4000  PRINT  USING  M12X.DD.2<4X.Z.4D>": J+Ja.R3< J-! > .R4(J-1 ) 

4010  OUTPUT  *®F  i  le3: J+Ja.R3< J-1 > ,R4< J-1 ) 

4020  NEXT  J 
4030  PRINT  '•  " 

4040  END  IF 

4050  IF  Nrun  MOD  5-0  AND  Mi -2  AND  Nrun<>30  AND  Iro-1  THEN 
4060  BEEP 

4070  PRINT  USING  "10X.""Set  porous-tube  flowmeter  reading  to  M",3D.D. 

. . :Set 

4080  END  IF 
4090! 

4100! 

4110! 

4120  IF  In-1  THEN 
4130  BEEP 

4140  INPUT  "DO  YOU  HAVE  ANY  MORE  DATA  <  1  -  YES  ,  0-N0>,?“  .Go_on 
4150  IF  Go_on=1  THEN  Repeat 
4160  ELSE 

4170  IF  Mi-2  AND  Nrun<30  THEN  Repeat 
4180  IF  Mi-1  AND  NrunOO  THEN  Repeat 
4190  END  IF 

4200  IF  Im- 1  THEN  PRINT  USING  "10X.DD,""  Data  runs  were  stored  in  fil 

Nrun . NewdataS 

4210  ASSIGN  KFilel  TO  * 

4220  ASSIGN  *File2  TO  - 

4230  ASSIGN  *File3  TO  * 

4240  END 
4250! 

4260!  THIS  SUROUTINE  CONVERTES  THERMOCOUPLE  VOLTAGE  INTO  TEMPERATURE 
4270! 

4280  SUB  Tvsv(T) 

4290  COM  /Ci/  C(7> 

4300  Sum-0. 

4310  FOR  1-0  TO  7 

4320  Sun-Sum+C< I ) »T *  I 

4330  NEXT  I 

4340  T-Sum 

4350  SUBEND 

4360! 

4370!  THIS  FUNCTION  CALCULATES  PRANDTL  NUMBER  OF  HATER  IN  THE 

4380!  RANGE  15  TO  45  DEG  C 

4390! 

4400  DEF  FNPrw(T) 

4410  Y-10  (1 .09976605-T*< 1 . 3749326E-2-T*<3 . 968875E-5-3. 4502SE-7*T> > > 

4420  RETURN  Y 
4430  FNEND 

4440! 

4450!  THIS  FUNCTION  CALCULATES  THERMAL  CONDUCTIVITY  OF  WATER 

4460!  IN  THE  RANGE  OF  15  TO  105  DEG  C 

4470! 

4480  DEF  FNKw<  T> 

4490  Y-.5625894  +  T*<2.296454SE-3-T*<! .509766E-5-4 . 0581 652E-3-T) ) 

4500  RETURN  Y 
4510  FNEND 

4520! 


"  PERCENT 


, 1 OA" : 
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4530! 

4540 ! 

4550! 

4560 

4570 

4580 

4530 

4600! 

4610! 

4620! 

4630! 

4640 

4650 

4660 

4670 

4680! 

n  y  _<  *j  . 

4700! 

4710 

4720 

4730 

4740 

05 

4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850! 
4860! 
4870! 
4880 
4890 
4900 
4910 
4920! 
4930! 
4940! 
4950! 
4960 
4970 
4980 
4990 
5000! 
5010! 
5020! 
5030! 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
5110 


THIS  FUNCTION  CALCULATES  SPECIFIC  HEAT  OF  WATER 
IN  THE  RANGE  15  TO  45  DEG  C 

DEF  FNCpw<T> 

Y-(4.21 ! 20358- T»< 2. 26826E-3-T«< 4. 4226 IE-5+2. 71428E-7*T>  >  >«I0G0 

RETURN  Y 

FNEND 

THIS  FUNCTION  CALCULATES  DENSITY  OF  WATER  IN  THE 
RANGE  15  TO  105  DEG  C 

DEF  FNRhow<T) 

Ro*999 . 52946+ T«< . 01 269-T»<5. 48251 3E-3-T+1 .234147E-5)) 

RETURN  Ro 
FNEND 

THIS  FUNCTION  APPLIES  CORRECTIONS  TO  THERMOCOUPLE  READINGS 

DEF  FNTemp<T.I) 

DIM  A(14) .B< 14) 

DATA  0.640533.0.573054.0.593101 .0.57298.0.56228.0.567384.0.569577 

DATA  0 . 553951 . 0 . 552008 . 0 . 566955 , 0 . 520998 . 0 . 522661 .0.531008.0.560783.0. 5524 

DATA  1 1 .8744.8.63163.9.39412,8.570246.3.299436.3.36677.3.04507.7.459766 
DATA  7 . 498928 . 7 . 9408 ,5.87072,5.391556,6.1 3399 . 6 . 48586 . 6 . 3262 2« 

READ  A< » )  .B<-> 

IF  I < 1 5  THEM 
T«T-<A< I)-3(I>+.001*T) 

ELSE 
T-T-.5 
END  IF 
RETURN  T 
FNEND 

THIS  FUNCTION  COMPUTES  THE  SPECIFIC  VOLUME  OF  STEAM 
DEF  FMVvst(T) 

V*58 . 4525588- T*< i .5 1508776-T  +  < . 0 1372746585-T-4 .253667 1 1E-5) > 

RETURN  V 
FNEND 

THIS  FUNCTION  CONVERTS  THE  VQLTAGE  READING  OF  THE  PRESSURE 
TRANSIDUCER  INTO  PRESSURE  IN  MM  HG 

DEF  FNPvsv< V) 

Y-1 . 1 103462+163. 36413*V 

RETURN  Y 

FNEND 

THIS  FUNCTION  CALCULATES  THE  SATURATION  TEMPERATURE  OF  STEAM  AS  A  FUNCTION 
OF  PRESSURE 

DEF  FNTvsp(P) 

IF  P<600  THEN 

T«31 .8776153+P»< .235854929-P*<3 . 66 1 3664E-4-P+2 . 41652372E-7) ) 

ELSE 

T-59 . 36562 +P»( .07379467-PM3. 15662E-5-P»6 .27246E-9) ) 

END  IF 
RETURN  T 
FNEND 


5!  20! 

5130!  THIS  -UNCTION  COMPUTES  THE  VISCOSITY  OF  WATER 
51  AO! 

5150  DEF  FNMuu(T) 

5  T  SO  Mu-1 .57S09473E-3-T*<3.5l 193576E-5-T*<3.5335816E-7-1 .3655861 15E-9-T ) > 
5170  RETURN  Mu 
5180  FNEND 
5190! 

5200!  THIS  FUNCTION  COMPUTES  THE  LATENT  HEAT  OF  VAPORIZATION 
5210! 

5220  DEF  FNHf g< T > 

5230  Hfg-2A97.7389-T*<2.2074+T«<1 .7079E-3-2 .8593E-6-T) > 

52A0  RETURN  Hfg 
5250  FNEND 
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B 


MODIFIED  WILSON-PLOT  METHOD 


1 000  . - . . . . . *•*»••»* 

1010!  FILE  NAME:  WILSON 
1020! 

1030!  THIS  PROGRAM  COMPUTES  THE  SIEDER-TATE 
1040!  COEFFICIENT  FOR  FLOW  IN  TUBES 

1050! . **-***-. . * . . * . . . * 

1060  COM  /Cc/  C<  7 ) 

1070  DIM  Vutl8>.Tci<!8>.Tco<!8>.Ts<18).Md<18>.Ta(18> 

1080  DATA  0 . 1 008609! .25727 .34389 .-767345.8295.78025595 .81 
1090  DATA  -9247486589 . 6 . 97688E M  ,-2.661 92E+ 1 3,3. 94078E  + 1 4 
1100  READ  C<*> 

1110  PRINTER  IS  701 
1i on  rptd 

1130  CLEAR  709 

1140  INPUT  "ENTER  MONTH.  DATE  AND  TIME  <HH:DD:HH:MM:SS>'\BS 
1150  OUTPUT  709 : "TD" : B$ 

1160  Di= . 01 4 1 
1170  Ai-PI»Di '2/4 
1130  Do-. 01 5875 
1190  L-.305 

1200  Km-21 .9 

1210  Ru- <Do-Di  )»Do/ (Km»< Do+-Di ) ) 

1220  Rf  i  =  0 . 

1230  Series:  ! 

1240  OUTPUT  709:"TD" 

1250  ENTER  709: AS 

1260  PRINT  USING  " 1 0X . ""Month .  date  and  time:  "".15A":AS 
1270  BEEP 

1280  INPUT  "ENTER  INITIAL  GUESS  FOR  SIEDER-TATE  COEFFICIENT". Ci 
1290  BEEP 

1300  INPUT  "ENTER  EXPGNENT  FOR  REYNOLDS  NUMBER". Xn 

1310  PRINT  USING  "1  OX . ""Ini t lal  guess  for  Sieder-Tate  coefficient  =  ”".Z.DD":Ci 
1320  PRINT  USING  "1  OX , ""Exponent  for  the  Reynolds  number  -  "”,Z.DD":Xn 

1330  BEEP 

1340  INPUT  "ENTER  THE  INPUT  MODE  < 1 -3054A ,2-FILE)" . Im 
1350  IF  Im- 1  THEN 
1360  BEEP 

1370  INPUT  "GIVE  A  NAME  FOR  THE  DATA  FILE" .D_f i leS 
1380  CREATE  BDAT  D_fileS.5 
1390  ELSE 
1400  BEEP 

1410  INPUT  "GIVE  THE  NAME  OF  THE  DATA  FILE" . D_f i leS 

1420  PRINT  USING  "1  OX , "”Fo 1  low  ing  analysis  uas  performed  for  data  stored  in  f  x I 
e  "".10A":D  files 
1430  BEEP 

1440  INPUT  "ENTER  THE  NUMBER  OF  DATA  RUNS  STORED". Mrun 
1450  END  I- 

1460  ASSIGN  jFile  TO  D_ftieS 
1470  BEEP 

1480  INPUT  "GIVE  A  NAME  FOR  PLOTTING  DATA  FILE", Dp  lots 
1490  CREATE  BDAT  OoiotS.5 
1500  ASSIGN  sFilep  TO  DplotS 
1510  BEEP 

1520  INPUT  "ENTER  ANALYSIS  TYPE  < 1 -HI .2-UI >" . It 

1530  IF  It-1  THEN  PRINT  USING  "1  OX . ""Analys is  type  -  Hi-METHOD . 

1540  IF  It-2  THEN  PRINT  USING  "10X, ""Analysis  type  -  Ui -METHOD""" 
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1550  K-0 

1 560  PRINT  "  " 

1570  PRINT  "  Oata  Vu  Tsat  Ti  To" 

1580  PRINT  "  *  <m/S)  (C)  <C>  <C>" 

1590  Repeat:! 

T600 ! 

1610!  RECORDS  THERMOCOUPLE  AND  PRESSURE  TRANSDUCER 
1620!  READINGS  AUTOMATICALLY  THROUGH  THE  HP  3054A 
1630!  AUTOMATIC  DATA  ACQUISITION/CONTROL  SYSTEM 
1 640 ! 

1650  IF  In- I  THEN 
1660  BEEP 

1670  INPUT  "ENTER  FLOWMETER  READING  <AS  A  PERCENT >” ,Fm 
1680  OUTPUT  709:"AR  AFO  ALO" 

1630  OUTPUT  722 : "F 1  R1  TI  Z1  FL1" 

1700  OUTPUT  709 : "AS  SA" 

•710!  READS  THERMOCOUPLE  FOR  WATFR  INLET 
1720  ENTER  722:  T< 0 ) 

1730  OUTPUT  709 : "AR  AF3  AL5" 

1740  OUTPUT  722 : "F 1  R1  TI  Z1  FL1" 

1750  FOR  1=1  TO  3 
1760  OUTPUT  709 : "AS  SA" 

1770!  READS  THREE  THERMOCOUPLES  FOR  WATER  OUTLET 
1780  ENTER  722;T<I> 

1790  NEXT  I 

1800  OUTPUT  709 : "AR  AF19  AL19" 

1810  OUTPUT  722: "FI  R1  TI  Z1  FL1" 

1820  OUTPUT  709: "AS  SA" 

1830!  READS  PRESSURE  TRANSDUCER 

1840  ENTER  722:P_voits 
1850  Pvap=FNPvsv<P_voI ts) 

1860  Ts<K)*FNTv«;p(Pvap) 

1870  Sum-0. 

1880  FOR  1=0  TO  3 

1890!  CONVERT  VOLTAGE  READINGS  TO  TEMPERATURE 
1900  CALL  Tvsv<  T< I) > 

1910  IF  1=0  THEN 

1920  Tci <K> -FN Temp <  T < 0 )  . 0) 

1930  ELSE 

1940  M-I  +  2 

1950!  APPLY  THERMOCOUPLE  CORRECTIONS 
1960  To=FNTenp( T< I) ,M) 

1970  Sum=Sum+To 

1980  END  IF 

1990  NEXT  I 

2000  Tco<K)-Sum/3. 

2010  ELSE 

2020  ENTER  ?Fi le: Ts<K> , Tci <K> . Tco(K) ,Fm 
2030  END  IF 

2040  Ta<K)-(Tci<K)+Tco(K>)*.5 

2050  Md<K> -66 . 36»Fm/< 100*60) 

2060!  COMPUTE  WATER-SIDE  VELOCITY 
2070  Vu(K)-Md<K)/(FNRhou(Ta(K) )-Ai> 

2080  BEEP 

2090  IF  Im-1  THEN  INPUT  "ARE  YOU  TAKING  MORE  DATA  < 1 -YES.O-NO>'>",Go_on 
2100  M-K  +  1 

2110  PRINT  USING  "!0X.DD.5<2X.DDD.DD)":M.Vw<K> .Ts(K) .Tci<K> ,Tco(K> 

2120  IF  Im-1  THEN  OUTPUT  «F i le ; Ts<K ) . Tc i <K> , Tco(K) ,Fm 

2130  K-K*1 

2140  IF  Im-1  THEN 


2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220! 
2230! 
2240 ! 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2530 

2590 

2600 

2610! 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 


IF  Go  on-1  THEN  Repeat 
ELSE  ' 

IF  M<Nrun  THEN  Repeat 

END  IF 

K-K-1 

J-0 

Jj-0 

PERFORM  ITERATION  TO  COMPUTE  SIEDER-TATE  COEFFICIENT 

Ssq*0 
Sx  =  0 . 

Sy  =  0  . 

Sxs-0 . 

Sxy =0 . 

J-J  +  1 
PRINT  “  - 

PRINT  ’*  Iteration  number  ■  ":J 

IF  J-1  OR  Jj-1  THEN  PRINT  ”  X  1/H] 

CFM 

FOR  1-0  TO  K 

0*Md<I)»<Tco(I)-TcitI))«4180. 

Lmtd*<Tco<I>-Tci<I)>  /L0G<  <Ts<I>-Tci(D>/(Ts(I>-Tco<I>)> 
Un-Q/<Lntd»PI-Do-L  > 

Unr-1 ./Un 

Rei  -FNRhowt  Ta<  I  > )  -Vw<  I  >  *Di/FNMu<  Ta(  I )  > 

Pru-FNPrw( Ta< I ) ) 

X-Re i  <-Xn)*Pru  (-.3333) 

IF  It-2  THEN 
Hi r-Unr-Di/Do 
Ku*FNKw<  Ta< I) > 

GOTO  2660 
END  IF 
IF  1-0  THEN 
Cf-1  . 

Ku*FNKui(  Ta(  I > ) 

Hi»i<w/Di*Ci-Rei'Xn*Pru'.3333*Cf 
Dt-Q/(PI-Di*L*Hi ) 

Cfc=<FNMu<  Ta< I) >/FNMu< Ta< I >+Dt >  > ‘ .  1  4  ’ 

IF  ABS<  <Cf c-Cf  >/Cfc)> .01  THEN 
Cf=<Cf+Cfc>/2. 

GOTO  2500 
END  IF 

Ho- 1 . /<Unr-Rw~Rf i*Do/D i-Do/<Di»Hi>) 

Hir-1/Hi 
Qo  =  0 
ELSE 

COMPUTE  1/HI 

Hi r- <  Unr- 1  /Ho*<  Q/Qo) ' .  3333-Rw-Rf  i-Do/Di )*Di /Do 

Dt-0/<PI*Di»L>*Hir 

Cf  *  <FNMu<  Ta( I > )/FNMu( Ta<I)+Dt))".!4 

END  IF 

IF  It-2  THEN 
CF-1  . 

H i c-Ku-C i /D i -Re i ~Xn*Pru' .3333*Cf 
Dt»Q/<PI*Di*L»Hic> 

Cfc-<  FNMu( Ta< I ) >/FNMu<  Ta< I ) +Dt >  >  ~ . J 4 
IF  ABS<<Cf-Cfc)/Cfc)>.01  THEN 
Cf-<Cf+Cfc>*.5 
GOTO  2680 


<TW-TB> 


2740  END  IF 
2750  END  IF 
2760  Qp-Q/<PI*Do*L) 

2770  X-X/Cf 

2780  IF  Jj-1  THEN 

2790  OUTPUT  ®FiIep;X,Hir 

2800  Hi re-A+B*X 

2810  Ssa-Ssq+(Hir-Hirc)'2 

2820  END  IF 

2830  IF  J-1  OR  Jj-1  THEN  PRINT  USING  "1  OX . Z . 7D , 3X , Z .7D . 4X ,DD .DD .2X . D .3DE .X . Z .5D 
"iX.Hir.Dt.Qp.Cf 

2840!  COMPUTE  COEFFICIENTS  FOR  LEAST-SQUARES  SCHEME 

•'850  Sx*=Sx+X 

2860  Sy*Sy+Hir 

2870  Sxs»Sxs+X~2 

2880  Sxy*$xy ►X«Hir 

?aqn  nfxt  t 

2900  N*K+ 1 

2910!  COMPUTE  THE  SLOPE  OF  THE  LEAST-SQUARES  LIME 
2920!  DEVELOPED  FOR  THE  WILSON  PLOT 
2930  B=(N*Sxy-Sy*Sx)/<N*Sxs-Sx*Sx> 

2940  Cic=*D  t  /  <  B*Kw ) 

2350  PRINT  USING  "! OX .""Intermediate  value  of  Sieder-Tate  coefficient  ■  "",Z.4D 
“  :Cic 

2960  IF  ABS<  <Ci-Cic)/Cic)>.0!  THEN 

2970  Ci-Cic 

2980  GOTO  2260 

2990  END  IF 

3000  IF  Jj=!  THEN 

3010  PRINT  USING  "10X .""Sieder-Tate  coefficient  »  "" , Z . 4D" :C ic 

3020  A-(Sy-B*Sx)/N 

3030  PRINT  USING  " 1 0X . ""Est imated  fouling  factor  »  "“ ,MZ .5DE . ""  <m  2-K/W) . :A 

3040  PRINT  "  Least-squares  line:" 

3050  PRINT  USING  " 1 3X . ""S lope  -  "".Z.5D":B 

3060  PRINT  USING  "1 3X ."" Intercept  -  "",MZ.5DE";A 

3070  ELSE 

3080  Jj-1 

3090  GOTO  2260 

3100  END  IF 

3110  PRINT  USING  "10X.""Sum  of  squares  *  ""  .D.5DE" :Ssq 
3120  ASSIGN  ■jFile  TO  - 
3130  ASSIGN  2'Fiiep  TO  » 

3140  IF  Im-1  THEN 
3150  BEEP 

3160  PRINT  USING  "10X.""NOTE:  "".DD.""  Data  runs  were  stored  in  file  "".SA";M,D 

_f i leS 

3170  END  IF 

3180  PRINT  USING  "10X.""N0TE:  "”.DD.""  X-Y  pairs  were  stored  in  file  "",10A":K+ 
1 .DplotS 
3190  BEEP 

3200  INPUT  "ARE  YOU  RUNNING  ANOTHER  SERIES  ( 1 >YES,0-NO>?".Go_on 
3210  IF  Go_on»1  THEN  Series 
3220  END 

3230  OEF  FMRhou(T) 

3240! 

3250!  THIS  FUNCTION  COMPUTES  THE  DENSITY  OF  WATER 
3260! 

3270  Ro*1006.35724-T«< . 774489-T**  <2 .262459E-2-T*3 . 03304E-4) ) 

3280  RETURN  Ro 
3290  FNEND 


X-Y  pairs  were  stored  in  file  "",J0A":K+ 


3300  DEF  FNPvsv(V) 

3310! 

3320!  THIS  FUNCTION  CONVERTS  THE  PRESSURE  TRANSDUCER  . 

3330!  READING  FROM  VOLTS  TO  PRESSURE  IN  MM  HG 
3340! 

3350  Y-1 . 1103482+163.3641 3»V 

3360  RETURN  Y 

3370  FMEND 

3380  DEF  FNTvsp<P) 

3390! 

3400!  THIS  FUNCTION  COMPUTES  THE  SATURATION  TEMPERATURE 

3410!  CORRESPONDING  TO  PRESSURE  IN  MM  HG 

3420! 

3430  IF  P<600  THEN 

3440  T-31 .8776158+PM .235854929-P*<3 . 66 1 3664E-4-P-2 . 41852372E-7) > 

3450  ELSE 

3460  T-59.36562+P»< .073794b/-P*<3. 15662E-5-P-S .2724GE-9) > 

3470  END  IF 
3480  RETURN  T 
3490  FNEND 
3500  DEF  FNTemp(T.I) 

3510! 

3520!  THIS  FUNCTION  APPLIES  THERMOCOUPLE  CORRECTIONS 
3530 ! 

3540  DIM  A<5) . 3<5) 

3550  DATA  0.640533.0.573054.0.593101.0.57298.0.567384.0.569577 
3560  DATA  1 1 .8744.8.63163.9.39412.8.570246.8.299436.8.36677 
3570  READ  A<  »)  .B<**> 

3580  T-T-(A<I)-B<I)*.001»T) 

3590  RETURN  T 
3600  FNEND 
3610  SUB  Tvsv<  T> 

3620  COM  /Cc/  C( 7) 

3630  Sum-O. 

3640  FOR  1-0  TO  7 

3650  Sam”Sum+C< I) *T ' I 

3660  NEXT  I 

3670  T-Sum 

3680  SUBEND 

3690! 

3700!  THIS  FUNCTION  COMPUTES  PRANDTL  NUMBER  FOR  WATER 
3710! 

3720  DEF  FNPr'i»<  T ) 

3730  Pr-10  < 1 . 099 76605 -T-*< . 0 1 3749326-T*<3 .968875E-5+3. 45026E-7-T) ) ) 

3740  RETURN  Pr 
3750  FNEND 
3760  DEF  FNMu<  T ) 

3770! 

3780!  THIS  FUNCTION  COMPUTES  THE  VISCOSITY  OF  WATER 
3790! 

3800  Mu- ! . 5037546575E-3-T*<  3 . 025732489E-5- T»(2 . S26439826E-7-T-8 . 1 860 1 937E- 101) 
3810  RETURN  Mu 
3820  FNEND 
3830  DEF  FNKu<  T) 

3840! 

3850!  THIS  FUNCTION  COMPUTES  THE  THERMAL  CONDUCTIVITY  OF  WATER 
3860! 

3870  Kwa-. 572183504477+1 ,52770121209E-3*T 
•3880  RETURN  Kwa 
3890  FNEND 
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PLOTTING  PROGRAM 


1000!  FILE  NAME :  PLOT 
1010  PRINTER  IS  705 
1020  BEEP 

1030  INPUT  "ENTER  MINIMUM  AND  MAXIMUM  X- VALUES" .Xnun.Xmax 
1040  BEEP 

1050  INPUT  "ENTER  MINIMUM  AND  MAXIMUM  Y-UALUES" . Ymin . Ymax 
1060  BEEP 

1070  INPUT  "ENTER  STEP  SIZE  FOR  X-AXIS" . Xstep 
1 080  BEEP 

10S0  INPUT  "ENTER  STEP  SIZE  FOR  Y-AXIS" . Ystep 
1100  BEEP 

1110  PRINT  "IN :SPt : IP  2300.1800.8300.6800:" 

1120  PRINT  "SC  0 . 1 00 .0.1 00 : TL  2,0:" 

I  1  ^  A  r>  r  *  r,  r\  t  t  \t  \ 

I  I  JU  Ol  A”  IUV  /  '  A I  lex  A  “  A(T|  1  1 1  / 

1140  Sfy*100/<  Ymax- Yru  n  > 

1150  PRINT  "PU  0.0  PD" 

1180  FOR  Xa=Xnin  TO  Xmax  STEP  Xstep 

1170  X=(Xa-Xmin>*Sf  x 

1180  PRINT  "PA" : X ."  . 0 :  XT:" 

1190  NEXT  Xa 

1200  PRINT  "PA  100.0-.PU:" 

1210  PRINT  "PU  PA  0.0  PQ" 

1220  FOR  Ya-Ymin  TO  Ymax  STEP  Ystep 
1230  Y-(Ya-Ymn)*Sfy 
1240  PRINT  "PA  0,":Y."YT" 

1250  NEXT  Ya 

1260  PRINT  “PA  0.100  TL  0  2" 

1270  FOR  Xa-Xnin  TO  Xmax  STEP  Xstep 
1280  X-<Xa-Xnin)*Sf x 
1290  PRINT  "PA":X,".1C O:  XT- 
1300  NEXT  Xa 

1310  PRINT  "PA  100,100  ?U  PA  100.0  PD" 

1320  FOR  Ya-Ymin  TO  Ymax  STEP  Ystep 

1330  Y-<Ya-Ynun>*Sfy 

1340  PRINT  "PD  PA  1 00 . " . Y . "YT" 

1350  NEXT  Ya 

1360  PRINT  "PA  100.100  PU" 

1370  PRINT  "PA  0.-2  SR  i .5.2" 

1380  FOR  Xa*Xmn  TO  Xmax  STEP  Xstep 
1390  X«<Xa-Xmin)*Sfx 
1400  PRINT  "PA" :X ." .0 : " 

1410  PRINT  "CP  -2.-1 :LB":Xa;"” 

1420  NEXT  Xa 

1430  PRINT  "PU  PA  0.0" 

1440  FOR  Ya*Ynin  TO  Ynax  STEP  Ystep 
1450  Y*< Ya-Ymin>*Sfy 

1460  PRINT  "PA  0 .  " Y  , 

1470  PRINT  "CP  -4 , - .25 :L8" ; va;"" 

1480  NEXT  Ya 
1490  BEEP 

1500  INPUT  "ENTER  X-LABEL".XlabelS 
1510  BEEP 

1520  INPUT  "ENTER  Y-LABEL".Y labels 

1530  PRINT  "SR  1.5.2:P'J  PA  50.-10  CP":-LEN<XlabelS)/2:"0:LB";XlabelS: 
1540  PRINT  "PA  -11.50  CP  0  .";-LEN<YlabelS>/2»5/6:"DI  0 , 1 :LB" : Y labe IS ; 
1550  PRINT  "CP  0.0  DI" 


1580 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
'660 
1670 
1680 
1690 
1700 
1710 
1720 
1 73C 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1330 
1840 
1850 
1360 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2030 
2090 
2100 
21  10 
2120 
2130 
2140 
2150 


Repeat : ! 

DEEP 

INPUT  "DO  YOU  WANT  TO  PLOT  DATA  PROM  A  FILE  < 1-YES. 0-N0>?",0k 

IF  01c -1  THEN 

BEEP 

INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE" .D_f i leS 

ASSIGN  -iFiIe  TO  D_fiieS 

BEEP 

INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER". Md 
BEEP 

INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIPS  STORED" .Npai rs 
BEE° 

INPUT  "SELECT  A  SYMBOL  FOR  THE  PLOTTER  < 1 . 2-+ . 3-c . 4-o .5-  > " ,Sy 
PRINT  "PU  DI" 

IF  Sy-1  THEN  PRINT  "SM*" 

IF  Sy-2  THEN  PRINT  "$M+" 

IF  Sy-3  THEN  PRINT  "$Mc" 

If-  <-  _  /,  T I  if"  I  <  DD  TklT  "CM*" 

r  oy  1  •*!»_•'«  rl\*ni  sj*  • 

IF  Sy *5  THEN  PRINT  "SM  " 

BEEP 

INPUT  "SELECT  MODE  (  I  *HN/H1  .2*HN< avg)/H1  )’* .  Jm 
IF  Md>Npaxrs  THEN 
FOR  1*1  TO  <Md-1> 

IF  Jm* 1  THEN  ENTER  *F i le : Xa. Ya. Yy 
IF  Jm-2  THEN  ENTER  sF i le ; Xa. Yy . Ya 
NEXT  I 
END  IF 

FOR  1-1  TO  Npairs 

IF  Jm-’  THEN  ENTER  sF i le ; Xa. Ya. Yy 

IF  Jn-2  THEN  ENTER  $F i le ; Xa. Yy . Ya 

X-<Xa-Xnin>*Sf x 

Y- <  Ya~ Ynin ) *Sf  y 

PRINT  "PA" .X. Y 

NEXT  I 

BEEP 

ASSIGN  ?File  10  « 

INPUT  "DO  YOU  HAVE  MORE  DATA  TO  BE  PLOTTED  < 1 -YES.O-NO)?" .Go_on 
IF  Go_on-1  THEN  Repeat 
END  IF 
BEEP 

INPUT  "DO  YOU  LIKE  TO  PLOT  THE  MUSSEL T  RELATION  ( 1  -  YES , 0*NO>?" .Go_on 
PRINT  "PU : SM" 

IF  Go_on-l  THEN 

FOR  Xa-Xmin  TO  Xmax  STEP  Xstep/50 
X-<Xa"Xmin>*$f  x 

IF  Jm* 1  AND  Xa>Xmin  THEN  Ya-Xa ' . 75- < Xa- 1 > ' . 75 
IF  Jm-2  AND  Xa>Xmin  THEN  Ya=Xa'<-.25> 

IF  Xa-Xmin  THEN  Ya-1 
Y«<Ya-Ymin>*Sfy 
PRINT  "PA" ,X . Y , "PD" 

NEXT  Xa 
BEEP 


INPUT  "MOVE  THE  PEN  TO  LABEL  THE  NUSSELT  LINE". Ok 
PRINT  "LBNusselt" 

END  IF 
BEEP 

INPUT  "DO  YOU  LIKE  TO  PLOT  KERN  RELATIONSHIP?". Yes 
IF  Yes-1  THEN 

FOR  Xa-Xmin  TO  Xmax  STEP  Xstep/20 


I 
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2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

233U 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2*40 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2530 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 


Ya-Xa  <-1/6> 

X*<Xa-Xnun)»Sf  x 
Y“(Ya-Ynun)»Sfy 
PRINT  "PA" .X . Y , ”PD" 

NEXT  Xa 
PRINT  "PU” 

BEEP 

INPUT  "MOVE  THE  PEN  TO  LABEL  KERN  RELATIONSHIP" .Ok 
PRINT  "LBKern :PU“ 

END  IF 

PRINT  "PU  PA  0.0" 

8FEP 

INPUT  "OK  TO  PLOT  CURVE  FOR  NEW  DESIGN  ( 1 -OK ,0-N0>?" .Ok 
IF  Ok • 1  THEN 

FOR  Xa-Xmin  TO  Xmax  STEP  Xstep/2 
IF  Xa* 1  THEN  Ya-1 

IF  Xa>!  THEN  Ya- < <Xa«5 ) ' . 75- <Xa*5- 1 >  ' . 75+5 > /6 
X*(Xa-Xmin)*Sf x 
Y*<  Ya-  Ymin)*Sfy 
PRINT  "PA" .X.Y."PD" 

NEXT  Xa 
PRINT  "PU" 

END  IF 
BEEP 

INPUT  "DO  YOU  LTKE  TO  PLOT  EISSENBERG  RELATION?” ,Go_on 
IF  Go  on* 1  THEN 

FOR  Xa-Xtnin  TO  Xmax  STEP  Xstep/10 
Ya-.6+.42*Xa  (-.25) 

X-CXa-Xmint-Sf  x 
Y*< Ya-Ymin)«Sfy 
PRINT  "PA"  ,  X  .  Y  ,."PD" 

NEXT  Xa 
PRINT  "PU" 

BEEP 

INPUT  "MOVE  THE  PEN  TO  LABEL  THE  EISSENBERG  LINE”. Ok 
PRINT  "LBE issenberg :PU" 

END  IF 
PRINT  "PU" 

BEEP 

INPUT  "DO  YOU  LIKE  TO  PLOT  THE  EXPTL  CURVE  <1-Y.0-N0>?".Go_on 

IF  Go  on* 1  THEN 

BEEP 

INPUT  "ENTER  THE  EXPONENT". Ex 
FOR  Xa*Xitun  TO  Xmax  STEP  Xstep/10 
Ya*Xa  < -Ex ) 

X*<Xa-Xrun)*Sf  x 
Y*<  Ya-Yrun)*Sfy 
PRINT  "PA" . X . Y , "PD". 

NEXT  Xa 
PRINT  "PU" 


BEEP 

INPUT 

PRINT 


PRINT  "PR  1 
PRINT  "LB": 


•MOVE  THE  PEN  TO  LABEL  THE  EXPTL  CURVE* 
•LBHN<avg>/H1-N" 

*PR  1.1" 

■i  R”  .  _r »  .  •••• 


END  IF 
BEEP 

INPUT  "DO  YOU  LIKE  TO  DRAW  A  STRAIGHT  LINE?".Go_on 

IF  Go  on*1  THEN 

BEEP 

INPUT  "ENTER  THE  SLOPE". SI 
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2770  BEEP 

2780  INPUT  "ENTER  THE  INTERCEPT" .Ac 

2790  FOR  Xa-Xnun  TO  Xmax  STEP  <Xmax-Xmm> 

2800  Ya“Ac+Si«Xa 

2810  Y-<Ya-Ynun)*Sfy 

2820  X-(Xa-Xnm)*Sfx 

2830  IF  Y<0  THEN 

2840  Xam-<  Ynun-Ac>/Sl 

2850  X*<Xam-Xnin)*Sf x 

2860  Y-0 

2870  END  IF 

2880  IF  Y> 1 00  THEN 

2390  Xam*  <  Ymax  Ac ; /S 1 

2900  X  =  (Xain-Xnir,;  »Sfx 

2910  Y  =•  1  00 

2920  END  IF 

2930  PRINT  "PA".X.Y."PD" 

2940  NEXT  Xa 

2950  END  IF 

2960  PRINT  "PU  SPO" 

2970  END 
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